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Genome-wide microarray analyses revealed that during biological activation of systemic 
acquired resistance (SAR) in Arabidopsis, the transcript levels of several hundred 
plant genes were consistently up- (SAR+ genes) or down-regulated (SAR - genes) in 
systemic, non-inoculated leaf tissue. This transcriptional reprogramming fully depended on 
the SAR regulator FLAVIN-DEPENDENT MONOOXYGENASE1 (FM01). Functional gene 
categorization showed that genes associated with salicylic acid (SA)-associated defenses, 
signal transduction, transport, and the secretory machinery are overrepresented in the 
group of SAR+ genes, and that the group of SAR~ genes is enriched in genes activated 
via the jasmonate (JA)/ethylene (ET)-defense pathway, as well as in genes associated with 
cell wall remodeling and biosynthesis of constitutively produced secondary metabolites. 
This suggests that SAR-induced plants reallocate part of their physiological activity 
from vegetative growth towards SA-related defense activation. Alignment of the SAR 
expression data with other microarray information allowed us to define three clusters of 
SAR+ genes. Cluster I consists of genes tightly regulated by SA. Cluster II genes can 
be expressed independently of SA, and this group is moderately enriched in H2O2- and 
abscisic acid (ABA)-responsive genes. The expression of the cluster III SAR+ genes is 
partly SA-dependent. We propose that SA-independent signaling events in early stages 
of SAR activation enable the biosynthesis of SA and thus initiate SA-dependent SAR 
signaling. Both SA-independent and SA-dependent events tightly co-operate to realize 
SAR. SAR+ genes function in the establishment of diverse resistance layers, in the direct 
execution of resistance against different (hemi-)biotrophic pathogen types, in suppression 
of the JA- and ABA-signaling pathways, in redox homeostasis, and in the containment 
of defense response activation. Our data further indicated that SAR-associated 
defense priming can be realized by partial pre-activation of particular defense 
pathways. 

Keywords: systemic acquired resistance, transcriptional profiling, salicylic acid, gene classification, gene 
regulation, defense priming 



INTRODUCTION 

Plants are equipped with a multi-layered immune system that 
employs constitutive and inducible defense strategies to antag- 
onize colonization by pathogenic microbes (Spoel and Dong, 
2012). In pathogen-inoculated plant tissue, conserved micro- 
bial structures (pathogen-associated molecular patterns, PAMPs) 
elicit PAMP- triggered immunity (PTI), a basal resistance response 
that contains the extent of infection by compatible pathogen iso- 
lates. Following recognition of race-specific pathogen effectors by 
plant immune receptors ("resistance proteins"), plants are able 
to activate effector-triggered immunity (ETI) that usually results 
in a hypersensitive response (HR) at inoculation sites and pro- 
vides effective local resistance to pathogens with a biotrophic or 
hemibiotrophic lifestyle (Jones and Dangl, 2006). Although ETI 
is associated with stronger local responses than PTI, the signal- 
ing networks underlying both resistance forms partially overlap 
(Tsuda et al, 2009). 



A localized microbial infection of a single or a few leaves 
can also immunize the rest of the foliage to subsequent infec- 
tion, a phenomenon known as systemic acquired resistance (SAR) 
(Fu and Dong, 2013; Shah and Zeier, 2013). Once activated, 
SAR provides enhanced resistance to a broad range of (hemi) 
biotrophic fungal, bacterial, and viral pathogens (Sticher et al., 
1997). The SAR response is initiated by microbes eliciting PTI or 
ETI at inoculation sites, and can also be triggered by localized 
leaf treatment with purified PAMPs. The mechanistic princi- 
ples leading to SAR induction by different types of bacterial 
pathogens and the resulting systemic immunization patterns are 
highly overlapping (Mishina and Zeier, 2007; Jing et al., 2011). 
For instance, compatible, PTI-inducing Pseudomonas syringae 
pv. maculicola ES4326 (Psm) and ETI-inducing Psm avrRpml 
elicit highly similar systemic responses in Arabidopsis thaliana 
(Arabidopsis) plants, including systemic accumulation of the SAR 
immune signals pipecolic acid (Pip) and salicylic acid (SA), and 
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enhanced systemic expression of a variety of classical SAR marker 
genes such as PATHOGENESIS-RELATED GENE 1 (PR1 ), PR2, 
and PR5 (Mishina and Zeier, 2006; Navarova et al, 2012). 

Effective long-distance communication between inoculated 
(1°) leaves and distant (2°, "systemic") leaves is required for 
the activation of SAR. Several plant-derived substances have 
been proposed to participate in SAR long-distance signaling 
(Dempsey and Klessig, 2012; Shah and Zeier, 2013). These involve 
the putative lipid transfer protein DEFECTIVE IN INDUCED 
RESISTANCE 1 (DIR1), the methyl ester of SA (MeSA), glycerol- 
3-phosphate (G3P), the diterpenoid dehydroabietinal, the dicar- 
boxylic acid azelaic acid, and the Lys catabolite Pip (Maldonado 
et al, 2002; Park et al., 2007; Jung et al, 2009; Chanda et al, 201 1; 
Chaturvedi et al., 2012; Navarova et al., 2012). The importance 
of several of these candidate signals for SAR induction in plants 
appears to depend on external parameters such as the light envi- 
ronment (Attaran et al, 2009; Liu et al., 2011; Navarova et al., 
2012). 

At the onset of SAR, the long-distance information released 
from 1° leaves is supposed to be perceived in 2° leaf tissue (Shah 
and Zeier, 2013), and a feedback amplification mechanism in 2° 
leaves that involves Pip and SA then ensures the activation of SAR 
(Navarova et al., 2012). The Lys aminotransferase AGD2-LIKE 
DEFENSE RESPONSE PROTEIN 1 (ALD1), whose expression is 
intensified in 1 and 2° leaves during SAR, is required for SAR 
activation (Song et al., 2004a,b). ALD1 generates the non-protein 
amino acid Pip, which has recently been identified as a critical 
metabolic SAR signal (Navarova et al., 2012). Pip-mediated resis- 
tance requires FLAVIN-DEPENDENT MONOOXYGENASE1 
(FMOl), another indispensable SAR module (Mishina and Zeier, 
2006; Navarova et al., 2012). Pip enhances both its own biosynthe- 
sis and downstream signaling in SAR via intensification of ALD1 
and FMOl expression, and systemic accumulation of the amino 
acid is required for the de novo synthesis of SA in 2° leaf tissue 
(Navarova et al., 2012). SA is a second critical SAR metabolite 
that is produced in plants from chorismate by ISOCHORISMATE 
SYNTHASE 1 (ICS1) (Wildermuth et al., 2001; Metraux, 2002). 
SA induces SAR-related gene expression via the downstream 
regulator NON-EXPRESSER OF PR GENES 1 (NPR1), a tran- 
scriptional co-activator and bona fide SA receptor (Wu et al., 
2012; Fu and Dong, 2013). SAR-induced plants therefore exhibit 
increased expression of a number of PR genes which is presumed 
to directly contribute to their state of increased disease resistance 
(Sticher et al., 1997). Additionally, SAR confers defense prim- 
ing, which enables plants to more effectively respond to future 
pathogen encounter (Jung et al., 2009; Navarova et al., 2012). 

The interaction between Arabidopsis and P. syringae repre- 
sents a useful model system to elucidate the molecular princi- 
ples underlying inducible plant immunity (Katagiri et al., 2002). 
Large scale transcriptional profiling in Arabidopsis has been 
used to better understand PTI- and ETI-associated defense net- 
works that are activated at sites of bacterial inoculation, and 
the mode of action of bacterial effectors to promote disease 
(Thilmony et al, 2006; Truman et al, 2006; Wang et al, 2008; 
Tsuda et al., 2009). Moreover, genome-wide microarray analy- 
ses have been employed to characterize Arabidopsis transcrip- 
tional responses induced by the synthetic resistance activator 



S-methyl-l,2,3-benzothiadiazole-7-carbothioate (BTH), which is 
often considered as a functional SA analogue (Wang et al., 
2006). A DNA microarray representing about 25-30 % of the 
Arabidopsis genes has also been used to monitor and analyse 
gene expression changes under different SAR-inducing conditions 
(Malecketal.,2000). 

In the current study, we aimed to characterize biologically- 
induced SAR in Arabidopsis at the whole genome level, classify 
SAR-regulated genes according to their function and regula- 
tion of expression, derive molecular and physiological charac- 
teristics of the SAR-induced state, and further clarify the role 
of FMOl in SAR. We therefore analysed the transcriptional 
changes that occur in upper 2° leaf tissue upon SAR induc- 
tion with Psm bacteria in lower 1° leaves by use of ATH1 
microarray chip analyses in the Arabidopsis wild-type and fmol 
mutant plants. These analyses revealed that SAR is associated 
with massive transcriptional reprogramming in systemic tis- 
sue that virtually fully depends on FMOl. Alignment of the 
SAR expression data with publicly available microarray infor- 
mation from defense-, stress-, and hormone-related experi- 
ments allowed us to obtain information about the regulation 
of genes that are up- and down-regulated during SAR. For 
instance, within the group of genes up-regulated during SAR, 
subgroups consisting of SA-independent, SA-dependent, and par- 
tially SA-dependent genes could be discriminated. Moreover, 
our evaluation indicated that overlapping and contrasting reg- 
ulatory principles exist for the induction of local resistance 
responses and SAR. Further, functional categorization of SAR- 
related genes suggested that, upon SAR induction, plants redi- 
rect some of their resources from vegetative growth towards 
defense. 

MATERIALS AND METHODS 

PLANT MATERIAL AND GROWTH CONDITIONS 

Arabidopsis thaliana (L.) Heynh. ecotype Col-0 and fmol mutant 
plants (T-DNA insertion line SALK_026163; Mishina and Zeier, 
2006) were grown in individual pots on an autoclaved mixture of 
soil (Klasmann, Beetpflanzensubstrat Typ R.H.P.16), vermiculite 
and sand (10:0.5:0.5) in a controlled environmental chamber (J- 
66LQ4, Percival, Boone, IA) within 9h day (9 AM to 6 PM; 
photon flux density 70 ixmol raT 2 s _1 ) and 15 h night periods in 
a relative humidity of 70%. Growth temperatures during the day 
and night period were set to 21°C and 18°C, respectively. SAR 
experiments were performed with 5-6 week-old plants exhibiting 
a uniform appearance. 

SAR EXPERIMENTS 

Overnight cultures of Pseudomonas syringae pv. maculicola 
ES4326 {Psm) were cultivated as described (Mishina and Zeier, 
2006). For SAR induction, plants were infiltrated between 10 AM 
and 11 AM into three lower (1°) leaves (typically leaf 7-9) with 
a suspension of Psm in 10 mM MgCi2 [optical density at 600 nm 
(OD 60 n) = 0.005]. Infiltration with 10 mM MgCl 2 served as the 
mock-control treatment. Upper (2°) leaves (typically leaf 10-12) 
were harvested and shock- frozen in liquid nitrogen two days after 
the treatment of 1° leaves. Together, three biologically indepen- 
dent SAR experiments were performed for microarray analyses 
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(see below), and two further biological replicates were performed 
for qPCR-based expression analyses of selected genes (Table 2). 

RNA ISOLATION 

Total RNA was isolated from frozen leaves using QIAzol® Lysis 
Reagent (Qiagen, http://www.qiagen.com/) following the manu- 
facturer's instructions. For each sample, two leaves from different 
plants that received the same treatment were used. 

QUANTITATIVE REAL-TIME PCR ANALYSIS 

RNA samples were reverse-transcribed with the Omniscript RT 
Kit (Qiagen) using 1 [ig of total RNA. 2.5 |xl of cDNA and 5 [l\ 
of SensiFAST™ SYBR No-ROX Kit (Bioline) were used in PCR 
reactions of 10 |xl total reaction volume containing 0.75 [iM 
gene-specific primers. The At4g26410 gene which is non- 
responsive to P. syringae-inoculation was used as a reference gene 
(Czechowski et al, 2005). The following primers were used to 
obtain the expression data for SAR experiments 4 and 5 (Table 2): 
PR 1 -FORWARD: 5'-GTGCTCTTGTTCTTCCCTCG-3', PR1- 
REVERSE: 5'-GCCTGGTTGTGAACCCTTAG-3', PR2-FW: 
5'-TCAAGGAAGGTTCAGGGATG-3', PR2-RV: 5'-GAGATTC 
ACGAGCAAGGGAG-3', PR5-FW: 5' - ATCGGGAGATTGC A 
AATACG-3', PR5-RV: 5'-ATGACCTTAAGCATGTCGGG-3', 
AGP5-FW: 5'-CTACTGAATCTCCACCAGCTC-3', AGP5-RV: 
5'-GAGGGAGACTCTGCTAACTG-3', CALM3-FW: 5'-GAC 
TGATGATAAATCGTTGGAG-3', CALM3-RV: 5'-CCCAACA 
AACTAAGCATCCT- 3' , LTPa-FW: 5'-GGTTCTACTTCTGACT 
CTCC-3', LTPa-RV: 5'-GTCCGTCTCCTTCTCCT-3', PBS3-FW: 
5'-TGCCTGCTCGAGTCGCAACC-3', PBS3-RV: 5'-TGGACTAA 
GCCACAGAGCAAATGGC-3', UGT76B1-FW: 5'-CTTTACA 
AGAGACTAAGGCAG-3', UGT76B1-RV: 5'-CACACCTATCT 
GTAACTTATCCC-3', 20GD1-FW: 5'-ACCAAATGCAGGTCA 
TAAGC-3', 20GD1-RV: 5'-TGAAGGGAAATAGAAAGTCGG- 
3', NIMIN-1-FW: 5'-AGTAAGAGAAGACGAAGAAGAG-3', 
NIMIN-1-RV: 5'-TCCGCCGTTAGATTTCCT-3', At4g26410- 
FW: 5'-GAGCTGAAGTGGCTTCCATGAC-3', At4g26410-W: 
5'-GGTCCGACATACCCATGACC-3'. The qPCR reaction was 
performed in triplicate in a Rotor-Gene Q apparatus (Qiagen) 
using the cycling program: 95°C for 2 min, followed by 40 
cycles at 95° C for 7 s, 60° C for 25 s, and finally 72° C for 3 min. 
The data was analyzed using the Rotor-Gene Q 2.0.2 software, 
setting the threshold of the normalized fluorescence to 0.1, which 
corresponded to the exponential phase of the fluorescence signal. 
The resulting Or and E values were used to calculate the relative 
mRNA abundance according to the AACj method (Livak and 
Schmittgen, 2001). The values were normalized to those of the 
reference gene and expressed relative to the mock-control sample. 

MICR0ARRAY ANALYSIS AND DATA EVALUATION 

For each SAR microarray experiment, RNA samples from at 
least 7 replicates for a particular condition (Col-0/mock, Col- 
0/SAR, fmol/mock, fmol/SAR) were mixed. The pooled RNA 
samples were quality- checked and expression profiling performed 
with the GeneChip® Arabidopsis ATH1 Genome Array using 
the 3' IVT express kit (Affymetrix) under accreditation condi- 
tions (DNAVision, Charleroi, Belgium). The quality of the used 
GeneChips was assessed and all the samples were hybridized, 



processed, and scanned in parallel, ensuring that samples could 
be directly compared to each other. Moreover, the raw microar- 
ray data was normalized using the RMA algorithm (Irizarry et al., 
2003a,b). 

Together, three biologically independent SAR microarray 
experiments were performed. Statistical analyses of the normal- 
ized expression values of the three biological replicates were 
performed using a two-sided Student's f-test. The large scale eval- 
uation of the microarray data was performed using Microsoft 
Excel® data sheets. Ratios of normalized expression values for 
Psm- and mock samples were calculated ([P/M]sar)> and genes 
were arranged and grouped according to the size of their averaged 
[P/M]sar values with the Excel® data sort function. Similarly, 
mean gene expression values gathered from other publicly avail- 
able microarray experiments were determined for each gene, 
and stimulus-to-mock ratios QS/MLtimulus) an d ratios of local 
expression values in Col-0 and mutant leaves following Psm inoc- 
ulation ([Col/mMtonf]p sm ) were calculated thereof (Table 6). The 
[P/M]sar ratios for each gene were aligned with the correspond- 
ing [S/MJstimuius and [Col/ mutant]p sm ratios using the "merge" 
function of FIRe, an Excel® macro designed for rapid microarray 
data analysis (Garcion et al., 2006). The detailed selection crite- 
ria for the categorization of genes are described in the main text. 
The Arabidopsis Information Resource (TAIR) functional cate- 
gorization tool was used to classify the genes according to Gene 
Ontology (GO) descriptions (http://www.arabidopsis.org/tools/ 
bulk/go/index.jsp). The SAR microarray data were deposited in 
the NASCArrays database (NASCARRAYS-703). 

DETERMINATION OF CAMALEXIN, ABSCISIC ACID. AND JASM0NIC 
ACID LEVELS 

For the time course analyses of the local and systemic levels of 
camalexin, JA, and ABA (Figure 8), three 1° leaves of Arabidopsis 
Col-0 plants were treated with Psm (OD 60 o = 0.005) or 10 mM 
MgCl2- At the indicated times after treatment, the treated (1°) 
leaves and three upper, non-treated (2°) leaves were separately 
harvested and frozen in liquid nitrogen. Individual samples 
consisted of 6 leaves from two plants. The determination of the 
leaf metabolite levels was performed by vapor-phase extraction 
and subsequent GC/MS analysis as described by Navarova et al. 
(2012). For quantification, metabolite peaks originating from 
selected ion chromatograms were integrated: camalexin (m/z 
200), ABA (m/z 190), and JA (m/z 224). The area of a substance 
peak was related to the peak area of dihydrojasmonic acid 
(m/z 156) or indole-3-propionic acid (m/z 130) for internal 
standardization. Experimentally determined correction factors 
were considered. 

RESULTS AND DISCUSSION 

FUNCTIONAL CHARACTERIZATION OF THE SAR STATE 

The bacterial phytopathogen Pseudomonas syringae pv. maculi- 
cola ES4326 (Psm) triggers a classical and robust SAR response 
in vegetatively growing Arabidopsis thaliana plants (Mishina and 
Zeier, 2006; Jing et al., 2011). SAR induced in Arabidopsis by 
compatible Psm or HR-inducing Psm avrRpml develops between 
day 1 and day 2 in the whole foliage after a localized leaf inoc- 
ulation has occurred, and the SAR response is fully established 
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2 days post inoculation (dpi) (Mishina et al., 2008). We aimed 
to broaden our understanding of the nature of the SAR state by 
gathering and analyzing the transcriptional changes that occur 
upon SAR establishment on the whole Arabidopsis genome level. 
Therefore, we inoculated lower (1°) rosette leaves of 5 week-old 
Arabidopsis vegetatively growing Col-0 plants with a suspen- 
sion of Psm (ODgoo = 0.005) and harvested upper, non-treated 
(2°) leaves two days after inoculation for RNA extraction. A 
mock-infiltration of 1° leaves with 10 mM MgCl2 served as a 
control treatment. Affymetrix expression analysis (GeneChip® 
Arabidopsis ATH1 Genome Array) was then used to compare 
the expression profiles in 2° leaves of mock- and Psm-treated 
plants. Considering the above-mentioned kinetics of SAR induc- 
tion in the Psm-Col-0 pathosystem, the selected time point at 
2 dpi allowed us to characterize the transcriptional reprogram- 
ming events when SAR has just fully established, but potentially 
not the earliest transcriptional events at the very onset of the 
response. We also included the fmo 1 mutant into the study which 
is fully compromised in SAR, systemic accumulation of SA, and 
systemic expression of characteristic SAR marker genes (Mishina 
and Zeier, 2006). 

In sum, we performed three independent SAR experiments 
that were conducted in the same growth chamber under identical 
light, temperature and humidity settings at different time periods 
with distinct batches of plants and pathogens. Each of the inde- 
pendent SAR experiments yielded one pooled RNA sample for 
both the mock-control and the SAR-induced state that was used 
for microarray analysis. The pooled RNA samples were derived 
from at least seven biological replicates within each SAR exper- 
iment and every biological replicate consisted of two 2° leaves 
from distinct plants. Thus, the gene expression samples result- 
ing from an individual SAR experiment exhibited a high intrinsic 
statistical validity. The 12 pooled RNA samples for microarray 
analysis (3 Col-0/mock, 3 Col-0/SAR, 3/moi/mock, 3/moi/SAR) 



were quality-checked and expression profiling performed using 
the 3' IVT express kit (Affymetrix) under accreditation conditions 
(DNAVision, Charleroi, Belgium). The raw ATH1 microarray 
data was normalized using the RMA algorithm and normalized 
expression values obtained (Irizarry et al., 2003a,b). 

To define genes systemically up-regulated during the SAR state 
in Col-0 plants (SAR + genes), we first calculated the ratios of 
the normalized expression values for the Psm- and mock-samples 
from individual SAR experiments for each gene [P/M]sar> deter- 
mined the mean values for the 3 ratios from different experiments 
for each gene, and selected those genes that were up-regulated by 
a factor of at least 3 on average. We further applied a two-sided 
f-test between the normalized expression values of the Psm- and 
those of the mock-samples for each gene and excluded those genes 
from our list with a P-value > 0.05. These two selection crite- 
ria yielded 305 genes out of the 22810 genes represented on the 
ATH1 GeneChip that were up-regulated by a factor of at least 3 
in a statistically significant manner among the 3 individual SAR 
experiments (Table 1A). The number of genes up-regulated upon 
SAR induction on average by a factor of 5 and 10 amounted to 
149 and 67, respectively. For the fmol mutant, not a single SAR+ 
gene existed based on these criteria (Table IB), corroborating our 
previous findings that functional FMOl is critical for the activa- 
tion of systemic defense responses and SAR (Mishina and Zeier, 
2006; Navarova et al., 2012). We also recognized that expres- 
sion of several genes was consistently suppressed following SAR 
establishment in Col-0 plants (SAR - genes), although the over- 
all degree of gene down-regulation was lower than the degree of 
up-regulation (Tables 1A,B). For instance, 17 and 276 genes were 
significantly down -regulated by a factor of at least 4 ( [P/M]sar < 
0.25) and 2 ([P/M]sar < 0.5), respectively. Again, not a single 
SAR~ gene was differently expressed in 2° leaves of Psm- and 
mock- treated fmol plants when following these selection criteria 
(Tables 1A,B). 



Table 1 | Number of genes up- (SAR + genes) and down-regulated (SAR genes) upon Psm-induced SAR in Arabidopsis Col-0 and fmol. 



A 


SAR+ (up-regulated) 


Col-0 


fmol 


B 


SAR (down-regulated) 


Col-0 


fmol 




[P/MIsar >10 


67 


0 




[P/M] SAR <0.1 


0 


0 




>5 


149 


0 




<0.25 


17 


0 




>3 


305 


0 




<0.5 


276 


0 


C 


SAR+ (up-regulated) 


Col-0 


fmol 


D 


SAR - (down-regulated) 


Col-0 


fmol 




[P/M] SAR >10 


145 


0 




[P/M] SAR <0.1 


2 


0 




>5 


295 


1 




<0.25 


50 


0 




>3 


547 


4 




>0.5 


700 


0 



(A,B) IP/M]sar symbolizes the mean value over SAR experiments 1, 2 and 3 of the ratios of the normalized expression values for Psm-samples divided by those of 
the mock samples. A two-sided t-test comparing the normalized expression values of the Psm- and the mock-samples was performed, and genes with P > 0.05 
were excluded. On the right of each table, the number of resulting (A) SAF& genes (for [P/M] S ar > 10, 5, or 3) and (B) SAR~ genes (IP/M1 S ar < 0. 7, 0.25 or 0.5) are 
given. (C,D) Final classification of SAR+ genes and SAPr genes. The data from the untypical SAR experiment 3 was excluded. (C) Only genes with mean IP/M1sar 
> 3 from SAR experiments 7 and 2 were considered for the SAR*' gene cluster. Genes whose individual IP/M1sar ratios in either experiment 7 or 2 were below 
2 were excluded. These selection criteria were taken as a basis for the final classification of genes into the SAR+ gene cluster (number of S4R+ genes in bold). 
Number of genes with mean IP/M] S ar > 5 and > 10 are also given. (D) Only genes with mean IP/MI S ar < 0.5 from SAR experiments 1 and 2 were considered 
for the SAFT gene cluster. Genes whose individual [P/MIsar ratios in either experiment 7 or 2 were above 0.67 were excluded (number of SAFT genes in bold). 
Number of genes with mean IP/M1 S ar < 0.25 and < 0.1 are also given. 
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We next examined the [P/M]sar expression ratios of indi- 
vidual genes from the distinct SAR experiments more closely. 
Remarkably, many genes that on average belonged to the most 
prominently SAR-induced genes were not up-regulated in a sta- 
tistically significant manner (P > 0.05) or only barely exhibited 
significantly increased expression values over all 3 independent 
experiments. Table 2 lists 10 representative genes from this group. 
Among them are genes such as PR1, PR2, PR5, and PBS3, which 
belong, according to previous results from other groups and our 
own findings (Sticher et al, 1997; Maldonado et al, 2002; Mishina 
and Zeier, 2006; Lee et al., 2007), to the most characteristic 
SAR + genes. We recognized that many of these genes were highly 
induced in SAR experiments 1 and 2, but showed only a mod- 
est or low degree of up-regulation in experiment 3 (Table 2). For 
instance, PR1 expression was induced by factors of 137, 183, and 7 
in SAR experiments 1, 2, and 3, respectively. Thus, albeit markedly 
up-regulated in each of the individual SAR experiments, the high 
quantitative differences between expression values of experiment 
3 compared to the two other experiments resulted in a P value 
larger than 0.05 for PR1. A similar trend was obvious for most 
of the other genes listed in Table 2. We consequently performed 
another two independent SAR experiments (experiments 4 and 5) 
and examined the expression characteristics of the 10 genes listed 
in Table 2 by quantitative real-time PCR analyses. Strikingly, in 
both newly conducted experiments, virtually all of the examined 
genes exhibited [P/M]sar ratio quantitatively similar to experi- 
ments 1 and 2. Furthermore, the [P/M]sar values derived from 
experiments 4 and 5 for the genes PR1, AGP5, UGT76B1, LTP- 
like, CALM3, and PBS3 were quantitatively much higher than the 
[P/M]sar values obtained from experiment 3. On the basis of this 
data and previous findings, we concluded that the bacterial inoc- 
ulation in experiment 3 only provoked a modest SAR response 
on the transcriptional level that is not representative for the SAR 
response that is generally observed. 

A subsequent systematic search for differently expressed genes 
in the control samples of experiments 1, 2, and 3 identified 
about 50 genes that differed in their normalized expression val- 
ues by a factor of 3 or more between experiment 3 and both 
other experiments. Interestingly, many of these genes represent 
central flavonoid pathway genes (Table 3). For instance, the nor- 
malized expression values of two main transcriptional regulators 
of anthocyanin biosynthesis, MYB90 and MYB75 (Borevitz et al., 
2000), were about two and one order of magnitude, respec- 
tively, higher for the experiment 3-samples than for the samples 
from experiment 1 or 2. Similar quantitative expression patterns 
existed for DFR and ANS, encoding two key enzymes of the 
anthocyanidin biosynthesis pathway, dihydrofiavonol reductase 
and anthocyanidin synthase, respectively, for the anthocyanin- 
5-O-glucosyltransferase gene UGT75C1, for the anthocyanin 
coumaroyltransferase gene A3GlcCouT and for the glutathione- 
s-transferase gene GSTF12 that is involved in anthocyanin accu- 
mulation (Saito et al., 2013). Therefore, the biosynthesis of 
anthocyanins apparently was activated in the plants employed for 
experiment 3 to a markedly higher extent than in the plants from 
the other two experiments. Anthocyanin production in leaves 
is a characteristic response of plants to unfavorable environ- 
mental conditions such as drought, nitrogen deficiency or high 
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Table 3 | Several central anthocyanin biosynthesis genes exhibit high normalized expression values in plants of SAR experiment 3 compared 
with experiments 1 and 2 (M1-3: mock treatment experiment 1-3; P1-3: Psm treatment experiment 1-3). 



Locus 


Name 


Gene description 






Expression value 








Ratio 








Ml 


M2 


M3 


P1 


P2 


P3 


M3/M1 


M3/M2 


At1g66390 


MYB90/PAP2 


PRODUCTION OF 
ANTHOCYANIN PIGMENT 2 


11.1 


17.3 


1281.1 


14.9 


12.5 


539.1 


115.4 


74.1 


At5g17220 


GSTF12 


Glutathione transferase, 
anthocyanin biosynthesis 


14.4 


19.6 


721.3 


15.1 


54.3 


284.1 


50.1 


36.7 


At5g42800 


DFR 


Dihydroflavonol reductase 


12.0 


12.7 


497.6 


12.1 


21.7 


117.9 


41.4 


39.2 


At4g22880 


ANS 


Anthocyanidin synthase, 
leucoanthocyanidin dioxygenase 


17.5 


31.4 


687.7 


17.8 


46.1 


134.4 


39.3 


21.9 


At4g32940 


y-VPE 


Vacuolar cysteine proteinase 


220.0 


315.7 


3362.8 


247.3 


315.2 


2297.1 


15.3 


10.7 


At5g54060 


UF3GT 


UDP-GLC-FLAVONOID 
3-O-GLUCOSYLTRANSFERASE 


15.2 


16.4 


192.1 


19.3 


18.3 


47.8 


12.7 


11.7 


At1g03495 


A3GlcCouT 


Anthocyanin 
coumaroyltransferase 


9.9 


15.4 


142.3 


10.4 


15.2 


45.6 


14.3 


9.2 


At4g 14090 


UGT75C1 


Anthocyanidin-5-O- 
glucosyltransferase 


10.6 


9.1 


86.7 


10.5 


9.2 


28.2 


8.1 


9.5 


At1g56650 


MYB75/PAP1 


PRODUCTION OF 
ANTHOCYANIN PIGMENT 1 


109.7 


117.4 


803.9 


28.7 


26.1 


481 .8 


7.3 


6.8 


At1g62710 


fr-VPE 


Vacuolar cysteine proteinase 


25.5 


36.4 


188.4 


17.0 


32.8 


104.5 


7.4 


5.2 


At3g51240 


F3'H 


Flavanone 3-hydroxylase 


155.6 


338.4 


1473.2 


128.6 


271.5 


544.6 


9.5 


4.4 


At5g 13930 


CHS 


Chalcone synthase 
(TRANSPARENT TESTA 4, TT4) 


26.8 


410.9 


1146.6 


93.6 


280.5 


334.9 


42.7 


2.8 


At1g65060 


4CL3 


4-coumarate-CoA ligase 


15.1 


31.9 


102.7 


22.0 


27.3 


22.0 


6.8 


3.2 


At5g05270 


CHI 


Chalcone isomerase 


10.8 


39.3 


84.3 


20.0 


23.4 


23.1 


7.8 


2.1 



light stress (Misyura et al., 2013). Although the plants used in 
experiment 3 did not exhibit a macroscopically obvious stress 
phenotype, it is likely that they had suffered an unexpected stress 
exposure prior to or in the course of the SAR experiment that 
resulted in the activation of anthocyanin biosynthesis. Mutual 
influences between pathogen defense signaling and abiotic stress 
pathways or the anthocyanin biosynthesis pathway do exist (Fan 
et al., 2009; Saijo et al., 2009), and it is possible that the moder- 
ate SAR reaction observed in experiment 3 was a consequence of 
such cross-talk events. 

Since the transcriptional data obtained from SAR experiment 3 
markedly differed from those resulting from the other four exper- 
iments, we decided to exclude this data set from further analyses. 
We defined new selection criteria to classify the SAR+ genes and 
only considered genes whose mean [P/M]sar ratios from experi- 
ments 1 and 2 were larger than 3, and whose individual [P/M]sar 
ratios in either experiment was at least 2. We thus ensured that a 
marked up-regulation of the selected SAR+ genes had taken place 
in each of the experiments 1 and 2. When following this proce- 
dure, all the characteristically SAR up-regulated genes listed in 
Table 2 fell into the category "SAR+ genes" which altogether con- 
sisted of 547 genes for the Col-0 wild-type (Table 1C). Thereof, 
295 and 145 SAR+ genes showed average [P/M]sar ratios larger 
than 5 and 10, respectively (Table 1C). To classify genes down- 
regulated during SAR, we selected genes with average [P/M]sar 
ratios lower than 0.5, whereby genes with rations higher than 0.67 
in either experiment 1 or 2 were excluded. This procedure yielded 
a group of 700 SAR~ genes, from which 50 genes had average 



[P/M]sar ratios lower than 0.25 and 2 genes exhibited average 
[P/M]sar ratios lower than 0.1 (Table ID). The new selection cri- 
teria pinpointed only 4 genes systemically up-regulated in fmol 
by a factor of at least 3 after Psm-infection. Moreover, not a single 
down-regulated gene was assigned for fmol (Tables 1C,D), again 
highlighting that fmol is virtually non-responsive to pathogen 
stimuli at the systemic plant level. 

As a first step to functionally characterize the two clusters of 
SAR + and SAR~ genes, we used the The Arabidopsis Information 
Resource (TAIR) functional categorization tool to classify the 
genes according to Gene Ontology (GO) descriptions which dis- 
criminates the three main classes "cellular component", "biologi- 
cal process" and "molecular function" (http://www.arabidopsis. 
org/tools/bulk/go/index.jsp). We thereby compared the SAR+ 
(547) and the SAR~ (700) gene cluster with all the genes rep- 
resented on the ATH1 Genechip (22810) (Tables 1C,D). In the 
class "cellular component", the functional categories "Golgi", 
"endoplasmatic reticulum (ER)", and "plasma membrane" were 
strongly overrepresented in the group of SAR+ genes. Moreover, 
the categories "cell wall", "extracellular", and "cytosol" were mod- 
erately overrepresented among the SAR+ genes, whereas the 
categories "plastid", "nucleus", and "ribosome" were underrep- 
resented (Figure 1 A). For the group of SAR~ genes, it became 
apparent that "chloroplast" "plastid", "cell wall", and "extracellu- 
lar" were highly represented categories (Figure 1A). These ten- 
dencies might indicate that cellular processes associated with the 
secretory apparatus and extracellular defenses are activated dur- 
ing SAR, and that certain activities occurring in chloroplasts and 
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chloroplast 
plastid 
mitochondria 
nucleus 
Golgi 
ER 

ri bosom e 
plasma membrane 
other membranes 
cell wall 
cytosol 
extracellular 
other cellular components 
other cytoplasmatic components 
other intracellular components 




ATH1 genes (N= 19696) 
SAR* genes (N = 492) 
SAR" genes (N = 646) 



B 

response to biotic or abiotic stimulus 
response to stress 
signal transduction 
transport 
transcription 
protein metabolism 
other metabolic process 
cellular organization and biogenesis 
developmental process 
DNA or RNA metabolism 
electron transport or energy pathways 
other biological processes 
other cellular processes 
unknown biological processes 




ATH1 genes (N = 19539) 
SAR + genes (N=521) 
SAR" genes (N = 664) 



kinase activity 
transferase activity 
transporter activity 
hydrolase activity 
transcription factor activity 
structural molecule activity 
other enzyme activity 
receptor binding or activity 
protein binding 
nucleotide binding 
DNA or RNA binding 
nucleic acid binding 
other binding 
other molecular functions 
unknown molecular functions 




ATH1 genes (N = 19075) 
SAR* genes (N=494) 
SAR genes (N = 61S) 




FIGURE 1 | Functional categorization of SAR+ and SAR" genes 
according to Gene Ontology (GO) descriptions 

(http://www.arabidopsis.org/tools/bulk/go/index.jsp). The 547 genes 
up-regulated during SAR by [P/M] SAR > 3 (SAR+ genes) (Table 1C), the 700 
genes down-regulated during SAR by [P/M]sar < 0.5 (SAR - genes) 
(Table ID), and all the genes (22810) represented on the ATH1 chip were set 



as input lists for the categorization algorithm. The depicted value on each of 
the x axes represents the quotient of "the number of genes annotated to 
terms of the respective categorization class" divided by "the total number of 
genes from the input list annotated to any term in this ontology (N)" in %. (A) 
Categorization class "cellular component'.' (B) Categorization class "biological 
process'.' (C) Categorization class "molecular function'.' 



plastids are reduced. Not unexpectedly, in the class "biological 
process", genes with GO annotations "response to (a)biotic stim- 
uli" and "response to stress" were typically found in the group 
of SAR up-regulated genes (Figure IB). Furthermore, the func- 
tional categories "signal transduction", "transport", and "protein 
metabolism" were highly represented in the group of SAR + genes, 
whereas the categories "developmental processes" and "DNA or 
RNA metabolism" were underrepresented. The category "electron 
transport and energy pathways" was relatively prominent among 
the SAR~ genes (Figure IB). Thus, the relative distributions of 
GO categories from the class "biological processes" indicated that 
stimulus- and stress-related signal transduction events, transport 
processes and protein metabolism are prominent features of the 



SAR state, whereas certain developmental, nucleic acid metabolic, 
and energy-related events might be reduced in SAR-induced 
plants. From the GO categories grouped into the class "molecular 
function", "kinase activity", "protein binding", "nucleotide bind- 
ing", "receptor binding", and "transporter activity" were overrep- 
resented in the SAR + gene cluster, whereas "DNA or RNA/nucleic 
acid binding" was poorly represented (Figure 1C). Again, this 
might emphasize the importance of signal transduction events 
such as protein phosphorylation and protein-ligand interactions 
for SAR. 

In the cluster of SAR~ genes, the relatively high abundance 
of the GO annotation "hydrolase activity" was one of the most 
obvious features (Figure 1C). Moreover, when examining the 
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specific functional annotation of genes in this cluster, two other 
trends became apparent. First, several genes belonging to the 
xyloglucan endotransglucosylase/hydrolase (XTH), the (fasciclin- 
like) arabinogalactan protein (AGP), the expansin-like protein, 
the extensin-like protein, and the polygalacturonase families 
were among the genes most strongly down-regulated during 
SAR (Table 4). Members of these gene families encode pro- 
teins associated with the extracellular matrix and/or the cell 
wall and have important functions in the rearrangement of 
cell wall components, wall loosening, cell elongation, and cell 
growth. For instance, members of the XTH family are involved in 
xyloglucan endotransglucosylation and/or in xyloglucan hydrol- 
ysis. These enzymatic activities can contribute to primary cell 
growth by restructuring and loosening the xyloglucan network, 
thereby enabling cell expansion (Rose et al., 2002). Two of the 
four XTH genes markedly down-regulated during SAR, XTH4 
and XTH31 (Table 4), code for proteins that have been exper- 
imentally identified as constituents of the cell wall proteome 
in Arabidopsis hypocotyls (Irshad et al, 2008). The strongly 
SAR down-regulated EXLA1 gene belongs to the expansin multi- 
gene family. The presence of EXLA1 protein in the Arabidopsis 
cell wall has also been experimentally verified (Irshad et al., 
2008). Expansins directly modify the mechanical properties 
of plant cell walls leading to turgor-driven cell extension (Li 
et al., 2002). Another class of extracellular proteins implicated 
in plant growth and development are the hydroxyproline-rich 
and highly glycosylated AGPs (Schultz et al., 2002). Several 
AGPs, among them the three fasciclin-like AGPs (FLAs) FLA8, 
FLA9, and FLA13, belong to the genes most highly down- 
regulated following SAR induction (Table 4). FLAs have, in addi- 
tion to predicted AGP-like glycosylated regions, putative cell 
adhesion domains known as fasciclin domains (Johnson et al., 
2003). Together, these examples indicate that the SAR state 
is associated with a marked down-regulation of various genes 
involved in cell wall modification, cell growth and development 
(Table 4). 



The cluster of SAR~ genes also contained several genes 
involved in glucosinolate and sinapoylester production (Table 5). 
These include genes encoding regulatory components or enzymes 
of indolic glucosinolate biosynthesis such as the MYB transcrip- 
tion factor MYB34 and the cytochrome P450 CYP79B3 that 
converts Trp to the indolic glucosinolate precursor indole-3- 
acetaldoxime (Glawischnig et al., 2004; Celenza et al., 2005), 
and of aliphatic glucosinolate biosynthesis such as the flavin- 
dependent monooxygenases FMOgs-oxi and FMOgs-ox3 that 
oxidize Met-derived methylthioalkyl glucosinolates to methyl- 
sulfinylalkyl glucosinolates (Li et al., 2008), or the 2-oxoglutarate- 
dependent dioxygenase AOP2 involved in the conversion of 
methylsulfinylalkyl to alkenyl glucosinolates (Table 5A) (Neal 
et al., 2010). Moreover, among the most prominently SAR down- 
regulated genes are SCPL8, SCPL10, and SCPL13 (Table 5B), 
encoding serine carboxypeptidase-like proteins that act as 
sinapoyltransferases to generate sinapoylmalate, sinapoylan- 
thocyanins, and 1,2-disinapoyl-glucose derivatives, respectively 
(Stehle et al., 2009; Fraser and Chappie, 2011). Sinapoylesters 
and glucosinolates are among the most abundant secondary 
metabolites produced in Arabidopsis in the course of nor- 
mal growth and development (Stehle et al., 2009; Sonderby 
et al., 2010). The reduced expression of genes involved in 
the constitutive production of major secondary metabolites 
in SAR-induced compared with control plants again supports 
the hypothesis that SAR represents a state of diminished 
vegetative growth. 

Together, the marked down-regulation of genes with pre- 
sumed roles in cell wall modification, cell growth and the con- 
stitutive production of secondary metabolites (Tables 4, 5), the 
overrepresentation of annotated chloroplast functions among 
SAR down-regulated genes (Figure 1A), and the strong up- 
regulation of stimulus-, stress- and defense-related genes during 
SAR (Figure IB) indicate that, compared with control plants, 
SAR-induced plants reallocate a part of their physiological activ- 
ity from vegetative growth towards particular defense-associated 



Table 4 | Genes associated with cell wall remodelling, cell extension, and growth belong to the most strongly down-regulated genes 
during SAR. 



Locus 


Name 


Gene description 


[P/MlsAR 


Rank among down- 
regulated genes 


At3g45970 


EXLA1 


Expansin-like A1 * 


0.14 


11 


At1g03870 


FLA9 


Fasciclin-like arabinogalactan-protein 9 


0.17 


13 


At2g45470 


FLA8 


Fasciclin-like arabinogalactan-protein 8 


0.18 


15 


At2g06850 


XTH4 


Xyloglucan endotransglucosylase/hydrolase* 


0.18 


16 


At3g06770 




polygalacturonase (pectinase)* 


0.20 


17 


At1g55330 


AGP21 


Arabinogalactan protein 21 


0.21 


22 


At2g 14890 


AGP9 


Arabinogalactan protein 9 /putative proline-rich protein 


0.21 


23 


At2g47930 


AGP26 


Arabinogalactan protein 26 


0.22 


26 


At1g12090 


ELP 


Extensin-like protein 


0.23 


32 


At5g44130 


FLA13 


Fasciclin-like arabinogalactan-protein 13 


0.24 


39 


At3g44990 


XTH 31 


Xyloglucan endotransglucosylase/hydrolase* 


0.24 


43 


At4g37800 


XTH7 


Xyloglucan endotransglucosylase/hydrolase 


0.25 


52 



Genes with asterisks encode proteins experimentally verified as cell wall constituents. 
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Table 5 | Genes involved in the biosynthesis of major constitutively produced secondary metabolites are down-regulated during SAR. 



Locus 



Name 



Gene description 



P/M Rank among down- 
regulated genes 



At5g60890 
At4g03060 



At4g03070 
At2g22330 



MYB34/ATR1 



AOP2 



At1 g65860 FMO GS-OX1 



At1 g62560 FMO GS-OX3 



AOP1 
CYP79B3 



MYB transcription factor, ALTERED TRYPTOPHAN REGULATION 1 , 0.25 51 
regulates indole glucosinolate biosynthesis 

2-oxoglutarate-dependent dioxygenase, ALKENYL HYDROXALKYL PRODUCING 2, 0.27 72 
aliphatic glucosinolate biosynthesis, conversion of methylsulfinylalkyl 
glucosinolates toalkenyl glucosinolates, not functional in Col-0 

FLAVIN-MONOOXYGENASE GLUCOSINOLATE S-OXYGENASE 1, aliphatic 0.34 163 
glucosinolate biosynthesis, conversion of methylthioalkyl glucosinolates to 
methylsulfinylalkyl glucosinolates 

FLAVIN-MONOOXYGENASE GLUCOSINOLATE S-OXYGENASE 3, aliphatic 0.36 215 
glucosinolate biosynthesis, conversion of methylthioalkyl glucosinolates to 
methylsulfinylalkyl glucosinolates 

2-oxoglutarate-dependent dioxygenase, similar to AOP2, 0.40 339 
possibly involved in aliphatic glucosinolate biosynthesis 

Cytochrome P450 monooxygenase, converts Trp to indole-3-acetaldoxime (lAOx), a 0.41 370 
precursor to IAA and indole glucosinolates 



At2g22980 SCPL13 

At2g23000 SCPL 10 
At2g22990 SCPL8/SNG1 



Serine carboxypeptidase-like (SCPL) protein; sinapoylglucose:sinapoylglucose 
sinapoyltransferase 

SCPL protein; anthocyanin sinapoyltransferase 

SCPL protein; sinapoylglucose:malate sinapoyltransferase 



0.11 



0.22 
0.30 



28 
100 



(A) Genes associated with glucosinolate biosynthesis. IB) Sinapoyltransferase genes involved in sinapoylester biosynthesis. 



processes that confer broad-spectrum disease resistance (see 
below). 

REGULATORY PRINCIPLES UNDERLYING SAR GENE EXPRESSION 

Having categorized defined groups of SAR up- and down- 
regulated genes, we next aimed to combine our SAR expression 
data with further transcriptional information to elucidate regula- 
tory principles that govern the SAR response. We used the FIRe 
software, an Excel® macro designed for rapid microarray data 
analysis (Garcion et al., 2006), to assemble the SAR expression 
data with expression data from other, publicly available microar- 
ray experiments describing the impact of various defense-, stress-, 
and hormone-related stimuli on gene expression in Arabidopsis 
plants. Information about the employed microarray data, the 
experimenters, and the experimental setup underlying each data 
set are summarized in Table 6. 

The information drawn from these microarray experiments 
is based on two distinct types of comparisons. In most exper- 
iments, Arabidopsis wild-type plants were exogenously treated 
with a chemical stimulus and the gene expression values of 
stimulus-treated plants or leaf tissue was compared to the val- 
ues resulting from an adequate mock-treatment. In this way, 
we could acquire information about the impact of exogenous 
application of defense and stress hormones or their derivatives 
[SA, methyl jasmonate (MeJA), abscisic acid (ABA)], oxidative 
stress (H2O2), the resistance-enhancing chemical BTH, which 
is often considered as an SA analogue (Lawton et al., 1996; 



Wang et al., 2006; Canet et al., 2010), and flg22-treatment, a 
22mer peptide corresponding to the elicitor-active domain of 
the bacterial PAMP flagellin (Gomez-Gomez et al., 1999), on 
gene expression (Table 6, Figures 2-7; microarray designations 
"BTH", "SA", "MeJA", "H 2 0 2 ", "ABA", "flg22"). Similarly, one 
experiment investigated the impact of Psm inoculation on gene 
expression (designation "Psm"). In contrast to our SAR microar- 
ray data (designation "SAR") that describes systemic changes in 
2°, non-inoculated leaves at 48 hours post inoculation (hpi), this 
experiment yielded information about the local changes in gene 
expression at the site of pathogen inoculation (1° leaves) at 24 
hpi. For all microarray experiments, we calculated the means of 
normalized expression values from the stimulus replicates and 
divided them by the means of the respective mock-values. This 
yielded, in analogy to the [P/M]sar ratios for the SAR experi- 
ment, stimulus to mock ratios [S/M] st ; mu i us that quantitatively 
indicated by which factors genes were differently expressed fol- 
lowing application of the exogenous stimulus compared with the 
mock-control in wild-type Col-0 plants (Figures 2-7). 

The microarray experiment 2 (Table 6) that investigated 
the impact of Psm leaf inoculation on local gene expression 
yielded two kinds of information: the ratio Psm/ mock in Col-0 
([S/M]p sm ; microarray 2a), and expression ratios of Psm-treated 
wild-type samples to different Psm-treated mutant samples 
([Col/mt/faHf]p sm ; microarray 2b). We evaluated the expression 
data from sid2 which is defective in ICS1 and, consequently, 
pathogen-induced SA production (Nawrath and Metraux, 1999; 
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Table 6 | Publicly available microarray data sets ("microarrays 1-7") used in this study. 



Microarray Name 



Experimenter 



NASCAR RAYS-703 Griebel, Attaran, 
Zeier 



2a 



2b 



NASCARRAYS-392 Wang, Dong 

(Wang et al., 2006) 



NASCARRAYS-454 Mitra, Glazebrook 
(Wang et al., 2008) 



NASCARRAYS-454 Mitra, Glazebrook 
(Wang et al., 2008) 



E-GEOD-3984 



Thibaud-Nissen 
(Thibaud-Nissen et al., 
2006) 



NASCARRAYS-174 Goda, Yoshida, 
Shimada 

(Goda etal., 2008) 



NASCAR RAYS-338 Mittler 

(Davletova et al., 2005) 

NASCARRAYS-176 Goda, Yoshida, 
Shimada 



(Goda et al., 2008) 



NASCARRAYS-123 Scheel, Brunner, 
Westphal 



Experiment description Designation 

Biological SAR, syringe infiltration of "SAR" 

lower leaves of 5 week-old Col-0 plants 

with Psm (OD 0.005), upper 

(non-treated) leaves harvested 2 d later, 

plants grown in soil under a 9/1 5-h 

light/dark cycle at 21/18 °C 

60>M BTH, spray-treatment of 4 "BTH" 
week-old Col-0 plants grown on soil 
under a 16/8-h light/dark cycle at 22° C, 
samples 24 h post treatment were 
considered 

Leaf inoculation (syringe infiltration) of "Psm" 
4-5 week-old Col-0 plants with Psm (OD 
0.002), inoculated leaves harvested 24 
hpi, plants grown in soil under a 12/12-h 
light/dark cycle at 22 °C 

Leaf inoculation (syringe infiltration) of "Col /mutant" 
4-5 week-old Col-0 or mutant plants with 
Psm (OD 0.002), inoculated leaves 
harvested 24 hpi, plants grown in soil 
under a 12/12-h light/dark cycle at 22 °C 

1 mM SA in 0.01 % Silwet, "SA" 
spray-treatment of 3-4 week-old, 
non-flowering Col-0 plants, leaf samples 
harvested 2 h post treatment 

7 day-old Col-0 seedlings grown in MS "JA" 
liquid medium under constant light at 
22°C were treated with 10p,M MeJA, 
leaf samples at 3 h post treatment were 
considered 

Application of 20 mM H2O2 to 5 day-old "H2O2" 

Col-0 seedlings grown on MS agar plates 

under constant light at 21-22° C 

7 day-old Col-0 seedlings grown in MS "ABA" 

liquid medium under constant light at 

22° C were treated with 10 [iM ABA, 

leaf samples at 3 h post treatment were 

considered 

Surface-treatment of leaves of 5 "flg22" 

week-old Col-0 plants with 1 mM flg22 

peptide, plants grown on soil at 22°C 

under a 8/16 hour light/dark regime, leaf 

samples 4 h post treatment were 

considered 



Depicted value 

Psm/mock = [P/M] SAR 



BTH/mock = [S/M] B 



Col-0-Psm/Col-0-mock = 
[S/M] Psm 



Col-0-Psm/ 'mutant-Psm 
= (Col/ mutant]p srn 



SA/mock = [S/MIsa 



JA/mock = [S/M] JA 



H 2 0 2 /mock = [S/M] H 202 



ABA/mock = [S/M] AB a 



flg22/mock= [S/M] flg22 



Sources, experimenters, relevant literature citations, and experimental descriptions are given. 
The designation of each experiment and the value depicted in the Figures 2-7 are also indicated. 



Wildermuth et al., 2001), SA insensitive nprl defective in the 
transcriptional co-activator and SA receptor NPR1 (Durrant and 
Dong, 2004; Wu et al, 2012), pad4 defective in the lipase-like 
defense regulator PHYTOALEXIN-DEFICIENT4 (Jirage et al, 
1999), ethylene-insensitive ein2 (Alonso et al., 1999), and JA 
insensitive coil defective in the JA receptor CORONATINE 
INSENSITIVE1 (Katsir et al, 2008). The [Col/tnutant] Psm ratios 



could be used to assess at the genetic level whether Psm-induced 
gene expression was dependent on SA accumulation (sid2), SA 
perception (nprl), JA perception (coil), and ET perception (ein2) 
(Figures 2-7). 

We now assembled the [P/M]sar ratios for each gene 
from the SAR experiment with the corresponding [S/M] or 
[Co\/mutant]p sm ratios from the other microarray experiments 
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microarray-experiment | 0 | 1 2a | 2b | : j 3 



2b 



gene cluster 


# 

genes 


SAR 
[PI**] 


BTH 

[S/M] 


SAR- genes 


524 


10.5 




SAR- genes 


670 






rest 


20575 





Psm 

[S/M] 


Psm 

[Col/ 
s/'d^] 


Psm 
[Col/ 
nprl] 


SA 

[S/M] 


Psm 
[Col/ 
pad4] 




.5 


3.0 


2.5 


2.5 


2.6 




.0 


0.7 


0.6 


0.8 


0.4 




.6 


1.2 


1.2 


1.0 


1.0 



Psm 

[Col/ 
ein2 



Psm 
[Col/ 
coil ] 



MeJA 

[S/M 



H 2 0 2 


ABA 


[S/M] 


[S/M] 



fig 22 

S/M] 





.9 1.5 1.0 3. 



.1 2.1 1.2 2. 



flg22 
[ S/M] 

1 



SAR genes 


# 

genes 


SAR 

[P/M] 


BTH 

[S/M] 


Psm 

[S/M] 


Psm 
[Col/ 
s;d2] 


Psm 
[Col/ 
npri] 


SA 

[S/M] 


Psm 
[Col/ 


Psm 

[Col/ 
e/'n2] 


Psm 
[Col/ 
coil ] 


MeJA 

[S/M] 


H 2 0 2 
[S/M] 


ABA 

[S/M] 


flg22 

[S/M] 


JA-activated 


76 


0.4 




6.6 


0.7 


0.5 


0.8 




1.1 


43.9 


KiUcl 


1.0 


1.5 


0.8 




























JA-repressed 


190 








0.6 


0.5 


0.8 




0.5 




0.7 


0.8 


1.3 


0.6 


JA-independent 


404 








0.8 


0.7 


0.8 




0.8 


0.9 


1.1 


0.9 


1.0 


0.6 




locally up 
24hpi 



systemically down 
48 hpi (SAR- genes) 
locally down 
24 hpi 



5 to 10 

2.5 to 5 



0.7 to 1.5 
0.4 to 0.7 
0.2 to 0.4 



FIGURE 2 | Mean values of [P/M] S ar, [S/M] stimulus , and [Col/ mutant]p sm 
ratios over groups of differently categorized genes. The numerical values 
are embedded in a heat map, and the legend on the right hand side depicts 
the value range assigned to each color of the heat map. The top row 
indicates the origin of the microarray data (Table 6). The selection criteria for 
the categorization of genes are detailed in the main text. (A) SAR+ genes, 
SAR - genes, and remaining ATH1 genes. Gene probes not unequivocally 



assignable to a single gene ("s_at"-probes) were removed so that the 
number of genes was slightly reduced in each group (Tables 1C,D). (B) 
SA-dependent (cluster I), SA-independent (cluster II), and remaining (cluster 
III) SAR+ genes. (C) JA-activated, JA-repressed, and JA-independent SAR~ 
genes. (D) Comparison of expression characteristics of SAR+ genes (row 1 ) 
and SAR - genes (row 3) with locally up-regulated genes ([S/M]p sm > 3, row 
2) and locally down-regulated genes ([S/Ml Psm < 0.5, row 4). 



using the "merge" macro of the FIRe program (Garcion et al., 
2006). This yielded an Excel® table in which the gene expres- 
sion information from all the experiments listed in Table 6 for 
the ATH1 genes was brought together. Some Affymetrix probes 
(labeled "s_at") hybridize to two or more related genes (Redman 
et al., 2004). These non-gene specific gene probes had been 
eliminated in some of the public microarray data used, and we 
consequently also deleted them from our merged Excel® list. The 
genes from our list were then grouped into different categories 
according to evaluation criteria outlined below (Figure 2), and 
excerpts of the Excel® data set relating to these categories are 
depicted in Figures 3-7. In addition to [P/M]sar ratios, the nor- 
malized expression values of the genes depicted in Figures 3-7 
for each of the three replicate SAR experiments is provided in an 
accompanying Excel® data file (Supplemental material). 



Before considering above-mentioned gene categories, how- 
ever, we discuss expression information of a specific example, 
the classical SAR marker gene PR1 (Figure 3, top row), to 
illustrate the gene regulatory information that we have drawn 
from the merged data set. PR1 is the most prominently up- 
rea ulated SAR+ gene ([P/M]sar = 160.2), and its expression is, 
as reported in previous studies (Delaney et al., 1995; Lawton et al., 
1996), enhanced by exogenous BTH ([S/M] BTH = 46.1) and SA 
([S/M]sa = 5.6) (Figure 3). Local inoculation with Psm also 
increases PR1 expression at 24 hpi, leading to an ([S/M]p sm ratio 
of 5.5 in the "Psm" experiment. The [Col/szd2]p sm (5.3) and the 
[Col/«pri]p sm (5.6) ratios reveal that PR1 expression upon Psm- 
inoculation in the two mutants does not exceed PR1 expression 
in Col-0 mock-control plants, indicating that P. syringae- induced 
PR1 expression fully depends on ICS 1 -mediated SA biosynthesis 
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Psm 


MeJA 


H 2 0 2 


ABA 


[Col/ 


[S/M] 


[S/M] 


[S/M] 


coil] 









flg22 
[S/M] 




gene description 



PATHOGENESIS-RELATED GENE 1 



NIM1-INTERACTING 1 

UDP-glycosy I transferase, 
conjugates isoleucic acid 

WRKY transcription factor 

acireductone dioxygenase (ARD) 
family 

GLUTAREDOXIN 13 

NUDIX HYDROLASE HOMOLOG 6, 
NADH pyrophosphatase activity 

WRKY transcription factor 

lysine/ histidine-like transporter 7 

P-loop containing nucleoside 
triphosphate hydrolase 
short chain 
dehydrogenase/ reductase 

WRKY transcription factor 

cupredoxin superfamily protein, 
putative L-ascor bate oxidase 
MILDEW RESISTANCE LOCUS O 2/ 
POWDERY MILDEW RESISTANT 2 

UDP-glycosyltransferase 

FLG22-INDUCED 
RECEPTOR-LIKE KINASE 1 
NONEXPRESSER OF PR GENES 1, 
SA signal transduction, SA receptor 



FIGURE 3 | Values of [P/M] S ar, IS/Mls,^^, and [Col/mufanf] Psnl 
ratios for individual genes of the group of SA-dependent (cluster I) 
SAR + genes (Figure 2B). The data for the NPR1 gene is also included 



into the list. The gene names of genes indispensable for SAR are 
depicted in red. The legend for the heat map representation is depicted 
in Figure 2 



and on NPR1 -mediated downstream signaling. Thus, exogenous 
SA is sufficient and endogenous SA accumulation following bac- 
terial inoculation is necessary to induce PR1 expression. By con- 
trast, induction of PR1 expression is independent of JA signaling, 
because exogenous MeJA ([S/M]ja = 1.2) does not elevate PR1 
levels, and the [Col/coz'i]p sm ratio equals 1.0, indicating iden- 
tical P. synttgae-induced expression of the gene in the Col-0 
wild-type and in JA-insensitive coil. Further, the [Col/pad4]p sm 
ratio equals 2.3, indicating an attenuated but not a fully compro- 
mised Psm-induced expression of PR1 in pad4 and thus a partial 
PAD4 dependency (Figure 3). Finally, PR1 expression occurs vir- 
tually independently of ET signaling ([Col/em2]p sm = 0.8) and 
is not stimulated by exogenous H2O2 ([S/M]h202 = 0.9), ABA 
([S/M]aba = 1.1), orflg22 ([S/M] flg22 = 1.1) (Figure 3). 

To draw information about regulatory principles of gene 
expression in the clusters of SAR up-regulated, SAR down- 
regulated, and remaining ATH1 genes, we first determined the 
mean values of [P/M]sar> [S/M] stimulus, an d [Col/ 'mutant]p sm 
ratios for all the genes from each category. Compared with the 
rest of the ATH1 genes, the SAR+ genes exhibited, in addi- 
tion to a strong average expression in 2° leaves of SAR-induced 
plants (mean [P/M]sar = 10.5), a marked average up-regulation 
in leaves of BTH-treated plants and in Psm-inoculated leaves 
(Figure 2 A). To a lesser extent, the average expression of these 
genes was stimulated by the SA pathway, H2O2, and flg22, and 



positively influenced by functional PAD4. Moreover, a small aver- 
age inducing stimulus of ABA on the expression of SAR + genes 
was obvious, and the JA- and ET- pathways had virtually no 
influence on the average SAR+ gene expression patterns. By con- 
trast, the genes down-regulated in SAR (mean [P/M]sar = 0.4) 
exhibited a completely different regulatory pattern. On average, 
these genes were markedly down-regulated by BTH and PAD4, 
and to lesser extent, by SA signaling, ET-signaling, and flg22- 
treatment. Remarkably, the average expression of the SAR~ genes 
was strongly stimulated by JA signaling ([Col/coz'i]p sm = 5.6) 
(Figure 2A). 

The positive effect of BTH, SA signaling, and PAD4 on the 
average expression of SAR+ genes is consistent with the facts that 
BTH induces plant resistance patterns similar to SAR (Lawton 
et al., 1996), that SA is a central signal for SAR (Wildermuth et al., 
2001; Mishina and Zeier, 2006), and that the PAD4 defense regu- 
lator is required for SAR establishment (Mishina and Zeier, 2006; 
Jing et al., 2011). One of the hallmarks of SAR is systemic SA 
accumulation at the onset of SAR (Metraux et al., 1990; Shulaev 
et al., 1995; Attaran et al., 2009), and increased levels of SA in 
2° leaves upon SAR induction is closely associated with increased 
expression of SAR- related genes (Shulaev et al., 1995; Mishina 
and Zeier, 2006). We therefore determined whether all or only 
a sub-fraction of the SAR+ genes are indeed up-regulated by 
SA. To categorize SA-regulated genes, we aimed at selecting only 
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gene description 



lipid-transfer protein/seed storage 2S 
albumin superfamily protein 



putative chitinase 



GH3 family protein, p-hydroxybenzoyl 
amino acid synthetase activity 



Toll-Interleukin-Resistance (TIR) 
domain family protein 



HEAT SHOCK PROTEIN 90-1 



GLUTATHIONE S-TRANSFERASE 22 



Calcium-binding EF-hand family protein 



BLUE-COPPER-BINDING PROTEIN, 
SENESCENCE ASSOCIATED GENE 14 



cytochrome P450 CYP71B15, 
camalexin biosynthesis 



cytochrome P450 family protein 



HEAT SHOCK PROTEIN 70/ HEAT SHOCK 
COGNATE 70-4 (HSC70-4) 



FLAVIN-DEPENDENT MONOOXYGENASE 1 



ATP binding cassette transporter 
ABCG40, abscisic acid transport 



ISOCHORISMATE SYNTHASE 1, 
salicylic acid biosynthesis 



WRKY transcription factor, 
regulates camalexin biosynthesis 



putative zinc finger transcription factor, 
RESPONSIVE TO HIGH LIGHT 41 



senescence-associated gene SAG13, 
short-chain alcohol dehydrogenase 



AGD2-like defense response protein 1, 
lysine aminotransferase, 
pipecolic acid biosynthesis 



cytochrome P450, indoleacetaldoxime 
dehydratase, camalexin biosynhtesis 



AvrRpt2-induced gene 



FIGURE 4 | Values of [P/M] SAR , [S/M] s , imu | us , and [Col/mutant] Psm ratios for individual genes of the group of SA-independent (cluster II) SAR 1 genes 
(Figure 2B). The gene names of genes indispensable for SAR are depicted in red. The legend for the heat map representation is depicted in Figure 2 



those genes whose induced local expression upon ftm-treatment 
was severely compromised in sid2. We therefore had to con- 
sider genes with increased [Col/sz'ii2]p sm ratios and first selected 
genes with [Col/sz'ii2]p sm ratios > 2. However, this criterion alone 
was not sufficient for selection because the genes strongly varied 
in their [S/M]p sm -ratios and thus their Psm-responsiveness. For 
instance, genes strongly Psm- up -regulated in Col-0 still exhibit 
considerable ftm-induced up-regulation in sid2 with the criterion 
[Col/ sid2]p sm ratios > 2 and thus are only weakly SA depen- 
dent. We consequently coupled the [Col/s7(i2]p sm ratio to the 
degree of ftm-responsiveness and defined that only genes with 
quotients of [Col/sid2]p sm I [S/M]p sm > 0.67 were taken. The 
combination of these two selection criteria provided a set of 
100 genes out of 541 SAR+ genes whose expression was locally 
Psm-inducible in a modest to strong manner and whose Psm- 
induced up-regulation was largely dependent on SID2/ICS1 and 
thus on endogenous SA. Moreover, most of these genes were 
also up-regulated by exogenous SA (Figures 2B, 3). To catego- 
rize SAR + genes independently expressed from SA, we assembled 
all the genes with a low [Co\/sid2]p sm ratio ([Col/W2]p snl < 1.5) 
that were up-regulated upon Psm inoculation by at least a fac- 
tor of 3. This selection yielded 156 SA-independent SAR+ genes 
(Figures 2B, 4). The remaining SAR+ genes (268) were grouped 



into a third category that mainly consisted of genes partly requir- 
ing SID2/ICS1 for Psm-induced expression (partly SA-dependent 
genes), or of genes not locally up-regulated by Psm at 24 hpi 
(Figures 2B, 5). Therefore, SAR + genes were categorized into 
three groups according to their SA-dependent expression: strictly 
SA-dependent genes (cluster I), SA-independent genes (clus- 
ter II), and cluster III genes predominantly consisting of genes 
with partial SA-dependency. Irrespective of their SA responsive- 
ness, the vast majority of SAR+ genes were also up-regulated 
by exogenous BTH, indicating that the action of the so-called 
"SA analogue" BTH on gene transcription is significantly broader 
than the action of SA itself (Figures 2B, 3, 4, 5). Similarly, genes 
down-regulated during biological SAR were generally down- 
regulated by exogenous BTH (Figures 2, 6). These tendencies 
indicate extensive overlap between the biologically-induced SAR 
state and the state of enhanced disease resistance after BTH appli- 
cation. Nevertheless, some differences between biological SAR 
and BTH-treatment existed for the transcription levels of indi- 
vidual genes not affected in biological SAR (Figure 2A, "rest"). 
Hereunder, 1.8 % of genes were positively ([S/M]bth > 3) and 
15.8 % negatively ([S/M] BTH < 2) regulated by BTH (specific 
examples are MAPKKK19, CYP94C1, UGT76E1, ACS2, DXL1, 
PR3; Figure 7). 
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gene description 



arabinogalactan protein 5 



Calcium-binding EF-hand family 
protein 



PATHOGENESIS- RELATED GENE 2, 
B-l,3-glucanase 



protein kinase superfamily protein 



PATHOGENESIS-RELATED GENE 5, 

osmotin / thaumatin-like protein 
C2 calcium/lipid-binding plant phos- 
phoribosyltransferase family protein 



WRKY transcription factor 



WRKY transcription factor 



PHYTOALEXIN-DEFICIENT 4, 
a/p-hydrolases superfamily protein 



glutaredoxin family protein 



ENHANCED DISEASE SUSCEPTIBILITY 
1, a/p-hydrolases superfamily protein 



WRKY transcription factor 



receptor-like protein kinase 1 



calmodulin binding protein 60-like.g 



PENETRATION!., 
syntaxin of plants 121 (SYP121) 



syntaxin of plants 122 



NPRl-like protein 3, SA receptor 



senescence-associated gene 101, 
acyl hydrolase 



mitogen-activated protein kinase 3 



calmodulin-binding protein-like 



FIGURE 5 | Values of [P/M]sar. [S/M] st j mu | us , and [Co\l mutant ]p sm ratios for individual genes of the group of partially SA-dependent (cluster III) SAR+ 
genes (Figure 2B). The gene names of genes indispensable for SAR are depicted in red. The legend for the heat map representation is depicted in Figure 2. 



On average, the [Col/ sid2]p sm ratios for the SA-dependent 
SAR+ genes, the SA- independent SAR+ genes, and the remaining 
genes amounted to 7.9, 1.1, and 1.6, respectively, reflecting strong, 
virtually absent, and moderate SA regulation of the respec- 
tive genes (Figure 2B). The [Col/«pri]p sm ratios (6.4, 1.4, 1.7) 
paralleled the [Coll sid2]p sm ratios, confirming above-mentioned 
influences of the SA pathway on the regulation of the differ- 
ent gene groups (Figure 2B). The [S/M]sa values exhibited a 
similar tendency as well (4.9 for SA-dependent, 2.1 for SA- 
independent, and 1.8. for remaining genes), although some genes 
of the SA-independent gene cluster showed a moderate respon- 
siveness to exogenous SA (Figure 2B). PAD3 is a typical example 
of a gene grouped into the SA-independent gene cluster which 
responded to exogenous SA (Figure 4). We reasoned, however, 
that in such cases, a [Col/ sid2]p sm ratio close to 1 would pro- 
vide a more meaningful criterion for SA-independency than an 
elevated [S/M]sa value, because the [Col/s;<i2]p sm ratio results 
from physiological differences in SA rather than from artificial SA 
differences caused by exogenous treatment. In general, however, 
the [Col/s!£?2]p sm values paralleled the [S/M]sa values remark- 
ably well: SA-dependent genes generally exhibited high values for 
both parameters (Figure 3), and the majority of SA-independent 
genes, as exemplified by CHI, FMOl, and SAG13, showed both 
[Col/sz'(i2]p snl and [S/M]sa values close to 1 (Figure 4). 



Noticeably, the average [P/M]sar ratios and the [S/M]fl g 22 
ratios were higher for the SA-dependent than for the SA- 
independent SAR+ genes, indicating a comparable high degree 
of up-regulation of SA-regulated genes upon both SAR induction 
and flg22-treatment (Figure 2B). When examining the expres- 
sion patterns of individual SA-dependent SAR+ genes, it became 
apparent that several genes such as FRK, LHT7, or UGT73D1 are 
strongly flg22-responsive whereas others are not at all (Figure 3). 
In addition, the average induction of SA-independent SAR + 
genes was higher than the induction of SA-dependent SAR+ 
genes in ftm-inoculated tissue, in H202-treated tissue, and in 
ABA-treated tissue (Figure 2B), indicating a more prominent 
stimulatory capacity of H2O2- and ABA-signaling on the expres- 
sion of individual members of the SA-independent compared 
with the SA-dependent gene cluster (Figures 2B, 3, 4). 

SAR+ GENES EXHIBITING TIGHT SA REGULATION (SAR+ GENE 
CLUSTER I) 

Representative examples of SAR+ genes tightly regulated by 
SA are PR1, the classical marker gene for SA-dependent 
defense gene activation (Nawrath and Metraux, 1999), NIMIN- 
J, UGT76B1, WRKY38, GRXS13, NUDX6, SDR3, WRKY70, and 
ML02 (Figure 3). NPR1, encoding a critical regulator of SAR 
that functions in SA perception and transcriptional activation of 
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gene description 



6.1 1.1 

0.8 0.7 

0.8 2.1 

2.5 0.6 

1.0 0.7 

0.9 1.4 

2.0 0.8 

1.0 0.9 

1.4 0.7 



JASMONATE-REGULATED GENE 21, 
2-oxoglutarate-Fe(II) oxygenase 



2-oxoglutarate-Fe(II) oxygenase 



benzoic acid/salicylic acid 
methyl transferase 



ornithine N-delta-acetyitransferase 



CORONATINE-INDUCED PROTEIN 1, 
CHLOROPHYLLASE 1 (CLH1) 



JASMONATE-ZIM-DOMAIN PROTEIN 5 
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isimilar to IAA amino acid conjugate hydrolase 
JASMONIC ACID RESPONSIVE 2, CORONATINE 
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defensin-like (DEFL) family protein 
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JASMONATE RESPONSIVE 1, 
JACALIN- RELATED LECTIN 35 (JAL35) 
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cytochrome P450, hydroperoxide lyase (non- 
functional in ecotype Col-0) 



GDSL-like lipase/acylhydrolase 
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JASMONATE-ZIM-DOMAIN PROTEIN 5 



POLYGALACTURONASE INHIBITING PROTEIN 

2 



MYB domain containing transcription factor, 
regulates anthocyanin biosynthesis 
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FIGURE 6 I Values of [P/M]sar. [S/M] s ,j mu | us , and [Col/ mutant]p sm ratios for individual genes of the group of JA-activated SAR genes (Figure 2C). 

The legend for the heat map representation is depicted in Figure 2. 



downstream genes, is only moderately up-regulated in 1° and 
2° leaf tissue upon inoculation ([P/M]sar = 2.6) and thus not 
assigned to the group of SAR+ genes. Nevertheless, its modest 
local up-regulation is also SA-dependent (Figure 3). 

Increased expression of PR1 upon SAR induction might 
directly contribute to resistance execution following fungal and 
oomycete pathogen attack, because PR1 proteins isolated from 
tobacco and tomato possess in vitro antifungal activity (Niderman 
et al., 1995). Moreover, overexpression of PR 1 in tobacco increases 
resistance to infection by the oomycetes Peronospora tabacina 
and Phytophthora parasitica var. nicotianae but has no pro- 
tective effect on tobacco mosaic virus or P. syringae infection 
(Linthorst et al., 1989; Alexander et al., 1993; Rayapuram et al., 
2008). Other SA-dependent SAR+ genes have a proven role 
in the activation of SA-associated defense responses that con- 
fers resistance to (hemi)biotrophic pathogens. For instance, the 
transcription factor WRKY70 has been recognized as a regula- 
tory node that positively regulates SA-related plant defenses and 



suppresses JA-mediated responses. Overexpression of WRBCY70 
increases basal resistance to P. syringae and to the powdery mildew 
Erysiphe cichoracearum, and results in the constitutive expression 
of SAR-related genes such as PR1, PR2, and PR5. Conversely, 
antisense suppression of WRKY70 or insertional inactivation 
leads to enhanced expression of JA-responsive genes and com- 
promises E. cichoracearum resistance (Li et al., 2004, 2006). In 
addition, WRKY70 acts in concert with WRKY53 and WRK46, 
two other SAR + genes that belong to the SA-dependent and 
the partial SA-dependent gene cluster, respectively (Figures 3, 5). 
This is reflected by the finding that wrky46/53/70 triple but not 
wrky70 single mutants exhibit attenuated basal resistance towards 
P. syringae (Hu et al., 2012). Another SA-regulated SAR compo- 
nent that positive regulates PR1 expression and is required for full 
basal resistance to P. syringae is the short chain dehydrogenase/ 
reductase SDR3. The metabolic function of SDR3 has not been 
elucidated yet (Hwang et al., 2012). Moreover, NUDX6, a mem- 
ber of the Nudix (nucleoside diphosphate linked to moiety X) 
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UOP-glucosyl transferase 76E12 



benzoic acid/salicylic acid methyl 
transferase 



MITOGEN-ACTIVATED PROTEIN KINASE 
KINASE KINASE 19 



1-AMINO-CYCLOPROPANE-l-CARBOX- 
YLATE SYNTHASE 2, ethylene biosynthesis 



1-DEOXY-D-XYLULOSE 5-PHOSPHATE 
SYNTHASE-LIKE 1 



TERPENE SYNTHASE 04, 
geranyllinalool synthase 



ALPHA-DIOXYGENASE 1, 
fatty acid oxygenation 
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NAC transcription factor family 



cytochrome P450 family protein; oxida- 
tion of 12-OH-JA-Ile to 12-COOH-JA-Ile 



cytochrome P450 family protein, desatu- 
ration of |}-sitosterol to stigmasterol 



PATHOGENESIS-RELATED GENE 3, 
basic chitinase, CHX-B 



UDP-glucosyltransferase 74E2, 
acts on indole-3-butyricacid 



cytochrome P450 family protein, bio- 
synthesis of the homoterpene TMTT 



UDP-glucosyltransferase 76E1 



cytochrome P450 family protein; 
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cytochrome P450 family protein; conver- 
sion of Trp to indole-3-acetaldoxime 



FIGURE 7 | Values of [P/M] SAR , [S/M] stimu | us , and [Col/ mutant] Psm ratios for individual genes strongly up-regulated at inoculation sites (24 hpi) but 
not in distal tissue (48 hpi) (Figure 2D). The legend for the heat map representation is depicted in Figure 2. 



hydrolase family that catalyze the hydrolysis of several nucleo- 
side diphosphate derivatives, not only acts in NADH metabolism 
in response to SA but also positively regulate SA-related defense 
responses (Ishikawa et al, 2010). MILDEW RESISTANCE LOCUS 
O 2 {ML02) belongs to a plant-specific family of genes cod- 
ing for membrane proteins that contain seven transmembrane 
domains. ML02 contributes to Arabidopsis resistance towards 
attack by the necrotrophic fungus Botrytis cinerea and partic- 
ipates together with other components to non-host resistance 
of Arabidopsis to the rice blast pathogen Magnaporthe oryzae 
(Humphry et al, 2010; Nakao et al., 2011). Interestingly, ML02 
has been recognized as target of the P. syringae type III effector 
HopZ2 which physically interacts with ML02. A mlo2 insertion 
line exhibits increased resistance to P. syringae, suggesting that 
the ML02/HopZ2-interaction is required for HopZ2-associated 
virulence (Lewis et al., 2012). 

NUDX6 and GRXS13 are two examples of SA-dependent 
SAR+ genes that appear to function in redox homeostasis dur- 
ing SAR. GRXS13 codes for a plant glutaredoxin which facilitates 
infection of Arabidopsis by B. cinerea (La Camera et al., 2011). 
Moreover, GRXS13 expression is critical to limit basal and pho- 
tooxidative stress-induced ROS production (Laporte et al, 2012). 
A redox-related function might also exist for At4g39830 which 
encodes a putative ascorbate oxidase (Yamamoto et al, 2005). 
Somewhat surprisingly, several SA-regulated SAR + components 
obviously reduce SA accumulation and/or SA signaling and there- 
fore appear to function in the containment of defense response 



activation during SAR establishment (Figure 3). For example, 
NIMIN-1 interacts with the SA receptor NPR1 in yeast-2-hybrid 
assays and functions as a negative regulator of SA-induced PR1 
expression (Weigel et al., 2005). The UDP- dependent glycosyl- 
transferase UGT76B1 can glycosylate the He catabolite isoleucic 
acid and thereby negatively influences SA accumulation (von 
Saint Paul et al, 2011). Beyond that, UGT76B1 exhibits in vitro 
glycosylating activity towards SA, and conversion of free, signal- 
ing active SA to glycosylated derivatives is supposed to attenuate 
SA signaling (Noutoshi et al., 2012). Finally, the transcription fac- 
tor WRKY38 negatively affects SA sensitivity and basal resistance 
to P. syringae (Kim et al., 2008). 

Together, SA-dependent SAR + genes can have distinct roles 
in the activation of defenses and resistance execution against 
different pathogen types (e.g. PR1, WRKY70, WRKY53, SDR3, 
NUDX6, ML02), down-regulation of the JA pathway (WRKY70), 
redox homeostasis (GRXS13, NUDX6), and containment of 
defense response activation after SAR establishment (NIMIN- 
1, UGT76B1, WRKY38) (Figure 9). The specific functions of 
several other SA-regulated SAR+ genes still remain to be 
clarified. 

SAR+ GENES THAT CAN BE EXPRESSED INDEPENDENTLY OF SA 
INCLUDE CRITICAL SAR ACTIVATORS (SAR+ GENE CLUSTER II) 

The group of SA-independent SAR+ genes contains at least 
three genes whose functions are necessary for SAR establish- 
ment: ALD1, FMOl, and ICS1 (Figures 4, 9). ALD1 encodes a Lys 
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aminotransferase (Song et al., 2004b) that mediates the biosyn- 
thesis of the Lys catabolite Pip, a critical SAR regulator (Navarova 
et al., 2012). Pip accumulation in 2° leaves of SAR-induced plants 
timely precedes SA accumulation, and Pip signaling requires the 
flavin-dependent monooxygenase FMOl for SAR induction. The 
function of Pip as a metabolic amplifier of defense signaling is 
crucial for the ICS 1 -mediated accumulation of SA in 2° leaves 
(Mishina and Zeier, 2006; Navarova et al, 2012). Therefore, at the 
onset of SAR, SA-independent signaling events obviously trigger 
the expression of the key SAR+ genes ALD1, FMOl, and ICS1 
that are required for systemic SA accumulation (Figure 9). ICS1 
can be regarded as a "bridge" between SA-independent and SA- 
dependent SAR signaling events because its up-regulation results 
in the de novo biosynthesis of SA (Wildermuth et al., 2001). It 
is important to note that a feedback amplification mechanism in 
2° leaves that involves ALD1, Pip, FMOl, ICS1, SA, and NPR1 
exists to ensure full SAR establishment (Navarova et al., 2012; 
Shah and Zeier, 2013). Therefore, in the context of SAR sig- 
naling and establishment in 2° leaf tissue, SA-independent and 
SA-dependent signaling process cannot be regarded as separately 
acting unities. For instance, although the SA-deficient sid2 mutant 
is able to accumulate wild-type like Pip levels in 1° leaves upon 
P. syringae inoculation, it accumulates markedly reduced Pip lev- 
els in 2° leaves (Navarova et al., 2012). Thus, Pip accumulation 
does occur independently from the capacity of SA biosynthesis at 
inoculation sites, but SA synthesis is required within the above- 
mentioned amplification cycle for the full accumulation of Pip in 
systemic tissue (Figure 9). Similarly, the systemic up-regulation of 
SA-independent genes might be substantially reduced in sid2, as 
has previously been shown for FMOl (Mishina and Zeier, 2006). 
Thus, the classification of SAR + genes as "SA-independent" is 
based on expression characteristics in 1° leaves (Table 6) and does 
not consider the necessity of SA production for SAR-associated 
signal amplification in 2° leaves. 

Another SAR+ gene up-regulated at inoculation sites in an SA- 
independent fashion is PBS3 (alias GH3.12, GDG1, and WIN3) 
(Figure 4). Null mutants of pbs3 exhibit severe defects in the 
induction of local resistance to bacterial infection and are signif- 
icantly but not fully compromised in SAR (Jagadeeswaran et al., 
2007; Lee et al, 2007; Nobuta et al, 2007). PBS3 acts upstream 
of SA in the induction of immune responses and encodes a GH3 
acyl adenylase that is able to conjugate 4-substituted benzoic 
acid derivatives to amino acids in vitro (Okrent et al., 2009). 
HSP90-1 and HSP70 are SA-independent SAR+ genes strongly 
up-regulated by H2O2 (Figure 4). HSP90-1 encodes a cytosolic 
isoform of the heat shock protein HSP90 that associates with 
the co-chaperones SGT1 and RAR1 to mediate the maturation 
of various nucleotide-binding domain and leucine-rich repeat 
containing (NLR)-type of resistance proteins. Gene knockouts 
of RAR1, SGT1 or HSP90 compromise resistance against vari- 
ous bacterial, fungal and viral pathogens (Kadota and Shirasu, 
2012). Arabidopsis SGTla but not SGTlb or SARI falls into 
the category of (SA-dependent) SAR + genes (Figure 9). HSP70 
(alias HSC70-4) represents one of four cytosolic HSC70 iso- 
forms that, similarly to HSP90, interact with the co-chaperone 
SGT1. Knockout of individual cytosolic HSC70 genes has no 
defense phenotype, but HSC70-1 overexpression compromises 



resistance to virulent and avirulent pathogens (Noel et al., 
2007). 

A typical local response of Arabidopis leaves to infection with 
necrotrophic or (hemi)biotrophic pathogens is the accumula- 
tion of the indolic phytoalexin camalexin (Glawischnig, 2007). In 
Psm-inoculated plants, camalexin is heavily produced in locally 
infected tissue but the phytoalexin does not accumulate sys- 
temically in 2° leaves (Figures 8A,B). However, SAR-induced 
plants are conditioned to more vigorously synthesize camalexin 
in response to subsequent pathogen encounter, and this priming 
effect is mediated by Pip (Navarova et al., 2012). The induction 
of camalexin biosynthesis at pathogen inoculation sites is associ- 
ated with a modest activation of genes coding for enzymes of the 
Trp biosynthetic pathway and strong induction of genes encoding 
cytochrome P450 monooxygenases involved in Trp catabolism 
towards camalexin: CYP79B2/3, CYP71A13, and CYP71B15 (alias 
PAD3) (Ren et al, 2008). In 2° leaves of SAR-induced plants, 
PAD3 expression is strongly and CYP71A13 expression is moder- 
ately up-regulated, and both genes belong to the SA-independent 
cluster (Figure 4). By contrast, Psm-inoculation does only lead to 
a local but not a systemic up-regulation of CYP79B2 (Figure 7), 
and CYP79B3 is neither locally nor systemically up-regulated. 
This expression pattern is consistent with the observed lack of 
systemic camalexin accumulation because CYP79B expression 
and thus the metabolic step from Trp to indole-3-acetaldoxime 
(IAOx) are not activated in 2° leaves of SAR-induced plants. 
However, enhanced systemic expression of CYP71A13 and PAD3 
pre-activates the pathway downstream of IAOx. This partial 
biosynthetic pathway activation can explain why camalexin accu- 
mulation is primed upon SAR induction and therefore can occur 
faster and more pronouncedly in challenged SAR-induced than in 
challenged control plants: key enzymes of the pathway are already 
expressed before the challenge inoculation takes place, and the 
ab initio expression of fewer components (for camalexin biosyn- 
thesis probably merely CYP79B2) is therefore required to provide 
the full enzymatic capacity for the biosynthesis of the metabolite 
(Figure 9). The WRK33 transcription factor controls the activa- 
tion of camalexin biosynthetic genes and camalexin production, 
and WRKY33 activity is regulated via a MAP kinase cascade 
involving MPK3 and MPK6 (Mao et al., 2011). WRKY33 and 
MPK3 but not MPK6 are up-regulated during SAR in a largely 
SA-independent manner (Figures 4, 5). 

As stated above, the average responsiveness to H2O2, ABA, 
and IA/ET is higher for SAR+ genes assigned to the SA- 
independent cluster than for those assigned to the SA-dependent 
cluster (Figure 2B). This is based on the fact that individ- 
ual SA-independent genes are strongly inducible by H2O2 (e.g. 
CYP81D8, WRKY33, BCB, HSP70, HSP90-1, ZAT12) or ABA (e.g. 
CHI, SAG13, GST22, HSP70), and that some genes are modestly 
inducible by JA/ET signaling (e.g. AIG2, GST22, PDR12, FMOl). 
Moreover, PDR12 encodes an ATP-binding cassette (ABC) trans- 
porter that mediates cellular uptake of ABA (Kang et al., 2010), 
and ZAT12 codes for a transcription factor implicated in ROS sig- 
naling (Davletova et al., 2005). This raises the question whether 
an early wave of ROS-, JA/ET-, and/or ABA-signaling would 
precede Pip- and SA-mediated establishment of SAR in 2° 
leaves (Figure 9). An earlier work on P. syringae-induced SAR 
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FIGURE 8 | Levels of camalexin, jasmonic acid (JA), and abscisic acid 
(ABA) in treated (1°) and non-treated distal (2°) leaves of Arabidopsis 
Col-0 plants inoculated with a suspension of P. syringae pv. 
maculicola [Psm; OD 0.005) or infiltrated with a 10 mM MgCI 2 solution 
(mock-treatment). Data represent the mean ± SD of at least three 
replicate samples. Asterisks denote statistically significant differences 



between Psm- and mock-samples (**P < 0.01; *P < 0.05; two-tailed f-test). 
(A,B) Camalexin levels at indicated times post treatment of 1° leaves in 
(A) 1° leaves and (B) 2° leaves. (CD) JA levels at indicated times post 
treatment of 1° leaves in (C) 1° leaves and (D) 2° leaves. (E) and (F) ABA 
levels at indicated times post treatment of 1° leaves in (E) 1° leaves and 
(F) 2° leaves. 



has reported the involvement of systemic ROS micro-bursts in 
the SAR regulatory network of 2° leaf tissue that requires an 
early oxidative burst in 1° inoculated leaves. These systemic 
microbursts were only observed in plants inoculated with avir- 
ulent but not with virulent P. syringae (Alvarez et al., 1998). 



However, the fact that both avirulent and virulent P. syringae are 
able to induce a typical SAR response argues against a critical role 
of microbursts in SAR activation (Mishina and Zeier, 2006, 2007; 
Attaran et al, 2009; Jing et al, 2011; Liu et al, 2011). Moreover, 
the occurrence of systemic microbursts during SAR has not been 
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FIGURE 9 | Events occurring in distal (2°) leaves of Arabidopsis plants in 
which SAR has been biologically activated by Psm inoculation in 1° 
leaves. SAR up-regulated (SAR+ genes) are depicted in red, SAR 
down-regulated genes (SAR - genes) depicted in green. Genes known to be 
indispensable for SAR activation are framed with a black line. The three 
clusters of SA-independent, SA-dependent, and partially SA-dependent SAR+ 
genes are illustrated. Processes leading to SAR establishment and functions 
of individual SAR-related genes or groups of genes are italicized. The 
alphabetical labels indicate a hypothetical order of events. (A) First, 
long-distance signals derived from inoculated leaves activate initial 
SA-independent signaling events. (B) Possible contributions of ROS, ABA, or 
JA/ET to these initial events are hypothetic or even doubtful. (C) A feedback 
amplification cycle (depicted as interconnected wheels) that requires the 



accumulation and the action of the two critical SAR metabolites pipecolic acid 
(Pip) and salicylic acid (SA) as well as the function of the 
flavin-monooxygenase FM01 establishes SAR. (D-F) Various events such as 
the activation of SA signaling, resistance induction, suppression of JA- and 
ABA-signaling, and partial pre-activation of camalexin biosynthesis occur. 
(G,H) SAR-induced plants are primed for early defense activation such as 
camalexin accumulation and defense gene expression. Small up arrows 
symbolize metabolite accumulation, small down arrows symbolize reduction 
of metabolite biosynthesis or of indicated physiological responses. Large 
arrows indicate the interconnection between the responses. Plus-signs 
symbolize activation, minus-sings repression. The indicated events are 
described in detail in the main text and summarized in the "Summary and 
conclusions" paragraph. 



confirmed by other studies. However, wounding and different 
abiotic stresses can trigger systemic ROS signaling that is depen- 
dent on the NADPH oxidase RBOHD (Miller et al, 2009). Miller 
et al., reported that the RBOHD-dependent signal can translocate 
from wounded to systemic tissue within minutes. This rapid sys- 
temic distribution of the ROS-related signal is not in accordance 
with the relatively slow establishment of SAR that essentially takes 
place between day 1 and day 2 after 1° pathogen inoculation for 
the Psm-Arabidopsis interaction (Mishina et al., 2008), and it is 
thus not clear whether ROS signaling indeed contributes to SAR 
establishment. 



Is the relatively high number of strongly ABA- and moderately 
JA/ET-responsive genes in the SA-independent group indicative 
for early ABA- and JA/ET-signaling events required for SAR 
establishment? Truman et al. (2007) have described an early JA 
response in 2° leaves of plants inoculated with very high titers 
[OD600 = 0.2] of avirulent P. syringae pv. tomato DC3000 (Pst) 
expressing the avrRpml avirulence gene (Pst avrRpml) that was 
associated with an increase of JA in petiole exudates collected 
from 1° leaves. These high initial titers of Pst avrRpml induce 
a rapid HR in leaves, and the necrotic disruption of leaf tis- 
sue goes hand in hand with the accumulation of JA and other 
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oxylipins (Mishina and Zeier, 2007; Zoeller et al., 2012). By con- 
trast, our experimental conditions include bacterial inoculations 
with lower densities (ODgoo = 0.005) of compatible Psm, which 
are not associated with an HR and elevations of JA in peti- 
ole exudates (Navarova et al., 2012). Although JA levels rise in 
ftm-inoculated leaves at about 48 hpi, systemic rises of JA are 
not detected between 16 to 48 hpi (Figures 8C,D). Considering 
these data and the above-mentioned 24 hpi to 48 hpi time win- 
dow in which SAR establishment takes place in ftm-inoculated 
Arabidopsis plants, it seems unlikely that JA signaling contributes 
to the up-regulation of the moderately JA/ET-responsive SAR+ 
genes of the SA-independent gene cluster. Irrespective of the 
existence of an early wave of systemic JA signaling after bacte- 
rial inoculation, its influence on SAR induction is supposedly 
negligible, because different Arabidopsis mutants compromised 
in JA biosynthesis or downstream signaling are SAR-competent 
(Attaran et al, 2009). 

The accumulation of ABA in Psm-inoculated leaves occurs 
more gradually and faster than the accumulation of JA 
(Figure 8E), and the kinetics of ABA production in leaves 
upon bacterial attack is similar for Psm- and Psf-inoculations 
(Figure 8E, Fan et al., 2009). Again, the accumulation of ABA, 
which negatively influences SA signaling and counteracts SAR 
(Fan et al., 2009), is confined to the site of bacterial inocu- 
lation and does not occur in 2° leaf tissue (Figure 8F). This 
metabolic data argues against a function for ABA as an initial trig- 
ger for the expression of SA-independent SAR + genes in 2° leaves 
during SAR. 

THE SAR+ GENE CLUSTER III CONTAINS GENES EXHIBITING PARTLY 
SA-DEPENDENT EXPRESSION 

The third cluster of SAR up-regulated genes contains genes 
that are essential for the initiation of the SAR process as 
well (Figures 5, 9). These include PAD4, ENHANCED DISEASE 
SUSCEPTIBILITY 1 (EDS1), SARD1, CBP60g, and NPR3 
(Mishina and Zeier, 2007; Truman et al, 2007; Zhang et al., 2010; 
Rietz et al., 2011; Fu et al., 2012). These genes are less tightly 
regulated by SA than the cluster I genes (Figure 3), because, 
unlike the latter, their [S/M]p sm ratios are (per definition) at 
least 1.5-fold higher than the corresponding [Col/ sid2]p sm val- 
ues (Figure 5). PAD4, EDS1, and another SAR+ gene product 
from cluster III, SAG101 (Figure 5), constitute a family of plant- 
specific hydrolase proteins that are critical regulatory components 
of plant basal resistance to (hemi)biotrophic pathogens and ETI 
triggered by a subset of resistance proteins (Wiermer et al., 2005). 
Apart from forming homodimers, EDS1 can interact with PAD4 
and SAG101, and the formation of different EDS1 homo- or 
heteromeric complexes is associated with distinct localization pat- 
terns in the cytoplasm and/or nucleus (Feys et al, 2005). Complex 
formation between EDS1 and PAD4 is required for the full estab- 
lishment of SAR (Rietz et al., 2011). Moreover, the EDS1-PAD4- 
SAG101 signaling complex also plays an important function in 
Arabidopsis post-invasion nonhost resistance to non-adapted 
powdery mildew fungi (Lipka et al., 2005). The [Col/pad4]p sm 
ratios illustrate that expression of the predominant fractions of 
SAR+ genes from all three distinguished clusters is more or 
less tightly regulated by PAD4, indicating an important function 



for PAD4 in overall SAR+ gene transcription (Figures 2B, 3, 4, 
5, 9). For instance, the expression of the critical SAR regulatory 
pathway genes ALD1, FMOl, and ICS1 are all partially depen- 
dent on PAD4 (Figure 5; Song et al., 2004a), and consequently, 
accumulation of the SAR regulatory metabolites Pip and SA are 
both positively regulated by PAD4 (Zhou et al., 1998; Navarova 
et al., 2012). Moreover, PAD4 tightly regulates the Psm-induced 
expression of two members of the plant-specific transcription 
factor gene family ACBP60, SAR-DEFICIENT1 (SARD1) and 
CALMODULIN -BINDING PROTEIN60G (CBPSOg), as well as 
the NPR1 homologue NPR3 (Figure 5). The SARD1 and CBP60g 
transcription factors are partly redundant in their function and 
activate pathogen-induced ICS1 transcription resulting in SA 
accumulation (Zhang et al., 2010; Wang et al., 2011). A sardl 
cbpg60g double mutant is therefore completely SAR defective 
(Zhang et al., 2010). NPR3 has been recently identified, besides 
NPR1 and NPR4, as a bona fide SA receptor (Fu et al., 2012; Wu 
et al, 2012). Therefore, the SAR+ gene cluster III contains central 
elements of both SA biosynthesis and SA downstream signaling. 

The SAR+ gene cluster III also contains the two WRKY 
transcription factor genes WRKY18 and WRKY40 (Figure 5). 
WRKY 18, WRKY40 and WRKY60 {WRKY60 is not a SAR+ 
gene) constitute a group of sequence-related WRKYs with com- 
plex and partly redundant roles in plant defense against different 
pathogen types. WRKY18 is required for biological SAR activa- 
tion and mediates a subset of NPRl-mediated responses (Wang 
et al., 2006). Whereas overexpression of WRKY18 alone increases 
resistance to P. syringae, simultaneous overexpression of WRKY 18 
and WRKY40 enhances susceptibility to the same pathogen (Xu 
et al, 2006). WRKY18/40/60 negatively regulate ABA signaling 
(Shang et al., 2010), and ABA signal transduction is also atten- 
uated by ARCK1, another cluster III SAR+ gene that encodes 
a receptor-like cytosolic protein kinase (Tanaka et al., 2012). 
Since ABA signaling can interfere with the SA pathway and thus 
attenuate plant defenses against (hemi)biotrophic pathogens (de 
Torres-Zabala et al., 2007; Fan et al., 2009), an impairment of 
the ABA pathway might ensure a robust SA response during SAR 
(Figure 9). Another cluster III SAR + gene is GRX480 encoding 
for a glutaredoxin that interacts with TGA transcription factors 
and negatively affects JA and ET signaling. This results in the sup- 
pression of expression of typically JA/ET-regulated genes such as 
PDF 1.2 (Ndamukong et al., 2007; Zander et al, 2012). 

The group III also contains the two SAR marker genes PR2 and 
PR5 which, similar to PR1, may function in the direct execution 
of disease resistance because their gene products exhibit antimi- 
crobial activity. Plant PR2 genes code for |3-l,3-glucanases, and 
purified |3-l,3-glucanases from pea have been shown to act syn- 
ergistically with chitinases in the degradation of fungal cell walls 
(Mauch et al., 1988). The PR5 protein family includes the basic 
osmotins whose members are homologous to the sweet-tasting 
protein thaumatin. Osmotin and other PR5 proteins have been 
shown to exhibit antifungal activity in vitro and in planta. For 
example, overexpression of tobacco osmotin in different plant 
species results in increased resistance to oomycete pathogens 
of the genus Phytophtora (Liu et al., 1994). Moreover, SAR+ 
group III contains PEN1, a critical determinant of Arabidopsis 
pre-invasion nonhost resistance to non-adapted powdery mildew 
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fungi. PEN1 codes for a plasma membrane-resident syntaxin 
which becomes recruited at sites of attempted fungal ingress and 
is implicated in a vesicle-associated resistance mechanism that 
prevents fungal penetration through epidermal cell walls (Collins 
et al., 2003; Bhat et al, 2005). The up-regulation of genes involved 
in non-host resistance, basal resistance to different pathogen 
types and ETI indicate that SAR simultaneously strengthens dif- 
ferent defense layers that make up the plant immune system 
(Thordal-Christensen, 2003). 

A SIGNIFICANT PORTION OF THE GENES DOWN-REGULATED DURING 
SAR EXHIBIT STRONG JA-RESPONSIVENESS 

A remarkable feature of the SAR~ genes is that their average 
expression is markedly activated via the JA signaling pathway 
(mean [Col/a>z'i]p sm = 5.6; Figure 2A). Taken the [Col/ani]p sm 
values as a basis, we divided the SAR~ genes into three categories: 
JA-activated genes ([Col/coi'i]p sm > 2; 76 genes), JA-repressed 
genes ([Col/cozi]p sm < 0.5; 190 genes), and JA- independent genes 
([CoVcoil] Psm > 0.5 and < 2; 404 genes) (Figure 2C). Although 
quantitatively the smallest group, the JA-activated genes most 
strongly influenced the average [Col/cozi]p sm ratio because most 
genes in this group have very high [Col/coz'i]p sm ratios (average 
[Col/coil ]p sm = 43.9). Apart from a few exceptions (e.g. PDF1.2), 
this is associated with high [S/M]/a ratios (average [S/M]ja = 
10.3) (Figures 2C, 6). Thus, both the [Col/coz'i]p sm and the 
[S/M]ja ratios indicate that the members of the JA-activated 
group of SAR~ genes are highly responsive to JA signaling. 

In fact, the JA-activated group of SAR~ genes consist 
of a series of genes typically regarded as marker genes of 
the JA pathway (Figure 6). Among them are JASMONATE- 
REGULATED 21 (JRG21) (Nickstadt et al, 2004), JASMONATE 
ZIM-DOMAIN PROTEIN 5 (JAZ5), JAZ9 (Thines et al, 2007), 
BENZOIC ACID/SA CARBOXYL METHYLTRANSFERASE 
1 (BSMT1) (Chen et al, 2003), N-ACETYLTRANSFERASE 
ACTIVITY 1 (NATA1) (Adio et al., 2011), CORONATINE- 
INDUCED PROTEIN 1 (COR1) (Benedetti et al., 1998), 
JASMONICACID RESPONSIVE 1 (JR1), JR2 (Leon et al., 1998), 
POLYGALACTURONASE INHIBITING PROTEIN 2 (PGIP2) 
(Schenk et al., 2003), ALLENE OXIDE SYNTHASE (AOS) 
(Laudert and Weiler, 1998), LIPOXYGENASE 2 (LOX2) (Bell 
and Mullet, 1993), and the plant defensin PDF1.2 (Ndamukong 
et al., 2007). This indicates that JA defense signaling is signifi- 
cantly reduced in the SAR-induced state. One of the key genes 
involved in JA biosynthesis, AOS (Laudert and Weiler, 1998), 
is markedly down-regulated during SAR (Figure 6), suggesting 
that the JA pathway could be attenuated already at the level of 
JA biosynthesis. However, since the experimentally determined 
levels of JA are similarly low in 2° leaves of mock-control and 
Psm-inoculated plants (Figure 8D), it is more likely that signaling 
events downstream of JA production are negatively affected 
during SAR. SAR is characterized by activated SA signaling 
(Figure 2A), and the well-established negative cross-talk between 
the SA- and JA-pathways might be responsible for the attenuation 
of JA responses (Spoel et al., 2003). As discussed above, molecular 
players such as the SA-activated transcription factor WRKY70 
and the glutaredoxin GRX480 could mediate the suppression of 
the JA pathway during SAR (Li et al, 2004, 2006; Ndamukong 



et al., 2007) (Figure 9). Moreover, the decreased expression of 
the SA methyl transferase BSMT1 that converts signaling active 
SA into inactive methyl salicylate (MeSA) supposedly counteracts 
deactivation of the SA signal (Attaran et al., 2009; Wu et al., 
2012) (Figure 9). SAR is associated with a reduced biosynthesis 
of constitutively produced metabolites such as glucosinolates or 
sinapoylmalates, as illustrated by the down-regulation of MYB34, 
a JA-inducible transcription factor that activates indolic glucosi- 
nolate production (Figures 6, 9) and by the decreased expression 
of several sinapoyltransferase genes (Table 5). Another group 
of secondary metabolites whose biosynthesis might be nega- 
tively affected upon SAR activation are anthocyanins, because 
MYB75 involved in the transcriptional regulation of anthocyanin 
biosynthesis (Borevitz et al., 2000) also belongs to the group of 
SAR- repressed and JA-activated genes (Figures 6, 9). It is likely 
that the production of inducible metabolites that negatively 
interfere with resistance to (hemi)biotrophic pathogens is also 
repressed during SAR. For instance, the JA-inducible acetyltrans- 
ferase NATA1 mediates the formation of N-8-acetylornithine 
from ornithine. Since natal knockout lines are more resistant to 
P. syringae than the wild-type, N-8-acetylornithine is supposed 
to negatively influence bacterial resistance (Adio et al., 2011). A 
reduced induction of NATA1 expression during SAR (Figure 6) 
could therefore counteract the presumed negative effect of N-8- 
acetylornithine on bacterial resistance. This would be consistent 
with the finding of Adio et al. (2011) that SA pre-treatment, 
which renders plants in an enhanced state of pathogen resistance, 
inhibits induced N-o-acetylornithine formation. 

SIMILARITIES AND DIFFERENCES OF DEFENSE ACTIVATION IN 1° AND 
2° LEAF TISSUE UPON BACTERIAL INOCULATION 

The signaling network underlying basal resistance to local infec- 
tion exhibits overlapping features to the signaling events that 
activate SAR because both forms of resistance share similar reg- 
ulatory factors such as ICS1, SA, NPR1, PAD4, EDS1, ALD1, 
Pip, and FMOl. ALD1 -mediated Pip production and FMOl- 
dependent transduction of Pip signaling do occur in both 
P. syn'ngae-inoculated and in systemic leaf tissue (Navarova 
et al., 2012). Pip functions as a mediator of defense ampli- 
fication in plants, and this fortification of defense responses 
is indispensable for the activation of SAR (Navarova et al., 
2012; Shah and Zeier, 2013). The Pip/FMOl -resistance path- 
way is also important for local resistance induction but the 
extent of its impact on basal resistance appears to vary with 
the attacking pathogen type (Song et al., 2004a; Bartsch et al., 
2006; Mishina and Zeier, 2006; Navarova et al, 2012). Further, 
SA accumulation and downstream signaling are common pro- 
cesses induced in 1° -inoculated leaf tissue and in distant 2°- 
leaves, and the activation of the SA pathway is required for 
both basal resistance to (hemi)biotrophic pathogens and SAR 
(Wildermuth et al, 2001; Mishina and Zeier, 2006; Spoel and 
Dong, 2012). The existence of common immune regulatory 
metabolites in 1° and 2° leaf tissue and the fact that most SAR+ 
genes are also up-regulated in inoculated tissue after pathogen 
encounter illustrates that overlapping signaling principles and 
defense mechanisms exist in inoculated 1° and in systemic 2° 
leaves (Figures 2A,B, 3, 4, 5). 
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Are there characteristic differences at the levels of defense 
metabolite production and gene activation in 1° and 2° leaf 
tissue? A first difference is of quantitative nature: Pip and SA 
accumulate to markedly higher levels in 1° than in 2° leaf tissue 
(Mishina et al, 2008; Navarova et al, 2012). In addition, the lev- 
els of a significant higher number of metabolites increase in 1° 
than in 2° leaf tissue after P. syringae inoculation (Ward et al., 
2010; Griebel and Zeier, 2010; Chanda et al., 2011; Navarova 
et al., 2012). As shown in Figure 8, the substances that strongly 
accumulate in 1° but not in 2° leaf tissue include JA, ABA, 
and camalexin. In addition, when comparing the number of 
genes up- or down-regulated in 1° and 2° leaves of ftm-treated 
plants in different experiments, it becomes obvious that the 
transcriptional changes at inoculation sites are much more pro- 
nounced than those in distant tissue (Figure 2D). From the 1921 
Arabidopsis genes assigned to be locally up-regulated at 24 h post 
ftm-inoculation (Figure 2D), 299 and 19 belong to the groups 
of SAR + and SAR~ genes, respectively. This implies that the 
expression levels of about 15.5 % of genes locally up-regulated 
at 24hpi do increase systemically after SAR establishment, but 
that the expression of the largest portion of locally induced genes 
essentially remains unchanged at the systemic level. Strikingly, the 
average [Co\/coil]p sm ratio in the group of locally up-regulated 
genes is high (6.4) compared to the same value for 2° leaves (1.1) 
(Figure 2D), indicating that JA signaling is strongly activated in 
1° but not in 2° leaf tissue, in which, as discussed in the previ- 
ous section (Figures 2C, 6), the expression of many JA- responsive 
genes is even reduced. The stimulation of JA signaling in 1° leaves 
can be triggered by the P. syringae phytotoxin coronatine, a struc- 
tural mimic of the signaling active JA derivative JA-Ile (Geng et al., 
2012), and by endogenously produced JA. In phases of the plant- 
bacterial interaction during which endogenous JA levels remain 
low (e.g. until 36 hpi for the ftrti-Arabidopsis inoculation experi- 
ment shown in Figure 8C), bacterial coronatine is presumably the 
major stimulus. For instance, Attaran et al. (2009) performed a 
comparative assessment of the formation of MeSA in Arabidopsis 
leaves induced by coronatine-producing and non-producing Pst. 
MeSA is generated by the BSMT1 -catalysed methylation of SA, 
and the BSMT1 gene is strongly JA- responsive (Figure 7). Until 
24 hpi, only the coronatine-producing but not the coronatine- 
deficient Pst strain elicited a significant formation of MeSA in 
inoculated plants, indicating that bacterial-derived coronatine 
rather than endogenous JA triggers the metabolic response in 
earlier phases of the interaction (Attaran et al, 2009). 

The strong activation of the JA pathway in inoculated leaves 
and the partial suppression of JA responses in distant leaves reflect 
a major difference between the hormonal status of 1° and 2° tis- 
sue, and this difference impacts the nature of defense responses 
in both tissue types. JA pathway activation at inoculation sites is 
at least partially causative for the more pronounced metabolite 
accumulation in 1° compared to 2° leaves because JA signaling 
induces the expression of a series of genes involved in metabo- 
lite biosynthesis. For example, the biosynthesis of many mono-, 
sesqui- and diterpenoids in plants is characteristically regulated 
via JA signaling (Arimura et al., 2008; Attaran et al., 2008). 
Indeed, among the most strongly up-regulated genes in Psm- 
inoculated tissue are JA-inducible genes involved in metabolism 



such as the SA methyl transferase BSMT1, the terpene synthase 
TPS4, several cytochrome P450 monooxygenases (e.g. CYP82G1, 
CYP94C1, CYP94B3, CYP79B2), and UDP-dependent glyco- 
syltransferases (UGT76E12, UGT76E1) (Figure 7). P. syringae- 
inoculated Arabidopsis leaves produce the Ci6-homoterpene 
(E,E)-4,8,12-trimethyl-l,3,7,ll-tridecatetraene (TMTT) in a 
TPS4-dependent manner (Attaran et al, 2008). TPS4 cataly- 
ses the first step of TMTT biosynthesis, i.e. the conversion of 
the diterpene precursor geranylgeranyl diphosphate to (E,E)- 
geranyllinalool (Herde et al., 2008). Subsequent formation of 
the Ci6-homoterpene from (E,E)-geranyllinalool by an oxida- 
tive cleavage reaction catalyzed by CYP82G1 completes TMTT 
biosynthesis (Lee et al., 2010). The two cytochrome P450 enzymes 
CYP94B3 and CYP94C1 are involved in the catabolic turnover 
of the signaling active jasmonate JA-Ile. CYP94B3 mediates the 
hydroxylation of JA-Ile to 12-hydroxy-JA-Ile and thereby inac- 
tivates hormone function (Koo et al., 2011). CYP94C1 then 
converts 12hydroxy- JA-Ile to the corresponding 12-carboxy- 
derivative (Heitz et al, 2012). Another strongly ftm-up-regulated 
gene whose expression is only moderately affected by the JA 
signaling pathway is 1-AMINO-CYCLOPROPANE-l-CARBOX- 
YLATE SYNTHASE 2 (ACS2), an ACS isoform involved in ethy- 
lene biosynthesis (Tsuchisaka et al., 2009). The activation of ET 
biosynthesis in P. syn'ngae-inoculated leaves is consistent with the 
microarray gene expression data, because ET-dependent genes are 
expressed more prominently in 1° (mean [Col/ez'«2]p sm = 1.7) 
than in 2° (mean [Col/em2]p sm = 1.1) leaves (Figures 2D, 7). 
Examples of genes partially regulated via ET signaling are 
ALPHA-DIOXYGENASE 1 (DOX1), encoding a fatty acid a- 
dioxygenases which converts linolenic acid and other fatty acids 
into their 2-hydroperoxy derivatives (Hamberg et al., 1999). 
Moreover, PR3 encodes a basic chitinase that is up-regulated in 1° 
inoculated but not in 2° tissue and regulated by JA/ET signaling 
(Figure 7; Zander et al., 2010). 

Another obvious difference between the transcriptional 
changes in 1° and 2° leaves following Psm inoculation is the 
stronger activation of ABA-responsive genes at inoculation sites 
compared to systemic tissue (mean [S/M]aba = 3.3 and 1.5 for 
1° and 2° leaves, respectively; Figures 2D, 7), which is consis- 
tent with the observation that ABA accumulation is restricted 
to bacterial inoculation sites (Figure 8E). The [S/M]aba values 
indicate that ABA signaling in 1° leaves contributes to the induc- 
tion of genes such as ANAC01 9, DOX1, CYP94B3, UGT74E2, and 
CYP710A1 (Figure 7). The transcription factor ANAC019 binds 
to a drought-responsive cis-element in the early responsive to 
dehydration stress 1 promoter, and overexpression of ANAC019 
in Arabidopsis provides increased drought tolerance (Tran et al., 
2004). This is one example of the fact that the local transcriptional 
responses following compatible P. syringae inoculation show large 
overlap with those occurring after drought or osmotic stress (de 
Torres-Zabala et al, 2007), presumably because strong bacterial 
proliferation in leaf tissue is associated with water deprivation 
and tissue necrosis. Activation of ABA signaling in infected tissue 
has also been interpreted as an active, effector-triggered virulence 
strategy of the pathogen, because ABA negatively interferes with 
SA signaling and therefore weakens one of the major pathways of 
plant defense to (hemi)biotrophic pathogens (de Torres-Zabala 
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et al, 2007; Fan et al, 2009). Moreover, the ABA- and H 2 0 2 - 
responsive UDP-depedent glucosyltransferase UGT74E2 that is 
also strongly up-regulated at inoculation sites has been implicated 
with the modulation of water stress responses. UGT74E2 glycosy- 
lates indole-3-butyric acid and therefore affects auxin homeosta- 
sis in plants (Tognetti et al., 2010), and modulated auxin signaling 
can result in disturbed plant immune responses (Truman et al., 
2010). The pathogen-induced CYP710A1 gene is also moderately 
ABA- and H2O2 -responsive (Figure 7), and CYP710A1 mediates 
the desaturation of the most common phytosterol in Arabidopsis, 
p 1 -sitosterol, to produce sitosterol (Morikawa et al., 2006), which 
strongly accumulates in P. syn'ngae-inoculated leaves, integrates 
into cell plasma membranes, and negatively affects plant defense 
and resistance to bacteria (Griebel and Zeier, 2010). 

SUMMARY AND CONCLUSIONS 

Based on gene expression data, metabolite data, and literature 
information, the present study aimed to contribute to a better 
understanding of the characteristics of the SAR-induced state in 
plants. Figure 9 summarizes main events occurring in 2° leaves of 
plants after biological SAR activation. These include: 

(1) The establishment of SAR in Arabidopsis in response to a 
localized leaf inoculation with the bacterial pathogen Psm is 
associated with a major transcriptional reprogramming in 
distant leaf tissue. Thereby, several hundred genes that are 
systemically up- (SAR + genes) and down-regulated (SAR - 
genes) can be distinguished. This extensive transcriptional 
reprogramming upon SAR induction is dependent on the 
SAR regulatory gene FMOl. 

(2) Functional categorization on the basis of GO annotations 
indicates that the SAR+ gene cluster is enriched in genes 
associated with stress responses, signal transduction, trans- 
port, and the cell secretory apparatus, whereas in the SAR~ 
gene cluster, genes associated with the chloroplast, cell wall 
loosening, cell extension, and the biosynthesis of constitu- 
tively formed secondary metabolites are over-represented. 
This suggests that, upon SAR induction, plants redirect 
some of their resources from vegetative growth towards 
defense-related processes. 

(3) Alignment of the SAR expression data with publicly avail- 
able microarray information has allowed us to classify the 
SAR-associated genes and analyse their expression charac- 
teristics. However, since the microarray data compared in 
our study originate from different laboratories, experimen- 
tal parameters such as plant age, growth conditions, the 
kind of treatment or the timing of sample collection var- 
ied between experiments (Tab. 6). For instance, whereas soil 
grown, 4 to 5 week-old Arabidopsis plants were used for 
the P. syringae-, the flg22-, the SA-, and BTH-treatments, 
5-7 day-old seedlings grown on MS medium were used 
for the JA-, ABA-, and H2O2 -treatments. Although these 
experimental differences might have impact on the expres- 
sion characteristics of individual genes, the predominant 
part of gene regulatory principles described in this work 
appears robust. This is exemplified by the findings that the 
JA- or SA-inducibility of genes was inferred from mutant 



analyses ([Col/ani]p sm or [Col/sid2]p sm ) and exogenous 
treatment ([S/M]ja or [S/M]sa) with considerable con- 
formity (Figures 3, 6). Moreover, the stimulus-dependent 
regulation of the (SAR-related) genes discussed in this study 
proved consistent with available literature data. 

(4) Based on the expression patterns of SAR-related genes in 
Psm-inoculated wild-type and sid2/icsl leaves and on the 
responsiveness of those genes to exogenous SA, we have 
categorized the group of SAR+ genes into clusters of SA- 
independent genes (cluster II), SA-dependent genes (cluster 
I), and partially SA-dependent genes (cluster III). 

(5) Albeit not congruent, extensive similarities of the transcrip- 
tional responses of Arabidopsis plants following biological 
SAR induction and treatment with the synthetic resistance 
activator BTH do exist. The so-called "SA analogue" BTH 
exhibits a broader effect on SAR-related gene expression 
than the endogenous defense signal SA. 

(6) The cluster of SA-independent SAR + genes contains the 
three critical SAR components ALD1, FMOl, and ICS1, 
which are indispensable for SAR establishment. ALD1 and 
FMOl are required for the biosynthesis and downstream 
signaling, respectively, of the immune regulator pipecolic 
acid, which mediates SAR activation via signal amplifica- 
tion (Navarova et al., 2012). ICS1 is involved in the de novo 
biosynthesis of SA (Wildermuth et al., 2001) and "connects" 
the SA-independent and SA-dependent phases of SAR. 

(7) In the initial stages of SAR establishment in 2° leaves, SA- 
independent signaling might precede and then activate an 
SA-dependent phase of SAR. Since both phases are required 
for the full activation of SAR, "SA-independent" and "SA- 
dependent" signaling events cannot be regarded as sepa- 
rately acting units but tightly co-operate to realize SAR. 
The average responsiveness to H2O2, ABA, and JA/ET is 
higher for SAR+ genes from the SA-independent than for 
those from the SA-dependent group. This might suggest a 
role of these stimuli in the early signaling stages of SAR in 
2° leaves. However, metabolite and mutant analyses rather 
argue against this possibility. 

(8) On average, SA-dependent SAR + genes exhibit a higher 
responsiveness to fig22 than SA-independent SAR+ genes. 
(Partially) SA-dependent SAR + genes function in the acti- 
vation or maintenance of distinct defense layers (non- 
host resistance, basal resistance, ETI), and in resistance 
execution against different pathogen types. This indi- 
cates that SAR heightens the plant immune system on 
several levels and illustrates the hallmark of SAR as a 
state of broad-spectrum resistance. Other SA-regulated 
SAR+ genes are involved in redox homeostasis and in 
the containment of defense response activation after SAR 
establishment. 

(9) Negative cross-talk between JA/ET- or ABA-signaling path- 
ways on the SA defense pathway is well-documented (Spoel 
et al, 2003; Fan et al., 2009). Several SAR+ genes are 
involved in the suppression of ABA- and JA/ET-signaling, 
suggesting that they can relieve inhibitory effects of these 
hormonal pathways on the SA pathway during SAR. 
Suppression of JA signaling during SAR also manifests itself 
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in the fact that many highly JA-responsive genes are among 
the most strongly down-regulated genes during SAR. 

(10) Overlapping defense principles exist for the induction of 
local resistance responses and SAR. However, the transcrip- 
tional and metabolic responses at sites of bacterial inocula- 
tion are generally more pronounced than those in systemic 
tissue. A major difference between the 1° inoculated and 
the 2° leaves relates to the stress hormonal status: whereas 
SA, Pip, JA, and ABA are produced at inoculation sites to 
largely high levels, only Pip and SA moderately accumu- 
late at the systemic level (Figure 8; Mishina et al., 2008; 
Navarova et al., 2012). Therefore, JA/ET- and ABA- triggered 
responses are strongly induced in 1° leaves, whereas these 
responses are not activated or even suppressed (see above) 
in the 2° leaves. 

(11) Pip accumulation during SAR primes plants to more quickly 
and effectively activate defense responses to subsequent 
pathogen encounter. A strongly primed defense response 
in Arabidopsis in SAR-induced plants is the Psm-triggered 
accumulation of the phytoalexin camalexin (Navarova et al., 
2012). SAR is associated with enhanced expression of the 
camalexin biosynthetic genes CYP71A13 and PAD3 but 
not CYP79B2, and therefore provides partial activation of 
camalexin production. This implicates that fewer compo- 
nents would have to be induced in a future pathogen chal- 
lenge to activate the whole response. Partial pre-activation 



of defense pathways might thus be a general mechanistic 
principle by which SAR-induced plants manage to accelerate 
defense responses when challenge-infected. 
(12) Varying environmental conditions can influence the magni- 
tude and particular mechanistic aspects of the SAR response 
(Shah and Zeier, 2013). The quantitative differences we have 
observed for the transcriptional SAR responses of exper- 
iment 3 compared with the responses in experiments 1 
and 2 were associated with higher overall leaf expression 
levels of major anthocyanin biosynthesis genes. Leaf antho- 
cyanin accumulation is a characteristic response to unfa- 
vorable environmental conditions (Misyura et al., 2013). 
An important task for future SAR research will be to 
systematically investigate to what extent and how other 
environmental issues and stress parameters influence SAR 
establishment. 
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